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. ABSTRACT 

High-resolution Hubble Space Telescope (HST) imaging observations of star cluster 
| systems provide a very interesting and useful alternative for stellar population anal- 

yses to spectroscopic studies with 8m-class telescopes. Here, we assess the system- 
atic uncertainties in (young) cluster age, mass, and - to a lesser extent - extinction 
and metallicity determinations, based on broad-band imaging observations with the 
HST. Our aim here is to intercompare the results obtained using a variety of com- 
monly used modelling techniques, specifically with respect to our own extensively 
tested multi-dimensional approach. Any significant differences among the resulting 
parameters are due to the details of the various, independently developed modelling 
techniques used, rather than to the stellar population models themselves. Despite the 
model uncertainties and the selection effects inherent to most methods used, we find 
that the peaks in the relative age and mass distributions of a given young (< 10 9 
yr) cluster system can be derived relatively robustly and consistently, to accuracies 
of a t = A(log(Age/yr)} < 0.35 and a M = A(log(M cl /M Q )) < 0.14, respectively, as- 
suming Gaussian distributions in cluster ages and masses for reasons of simplicity. 
The peaks in the relative mass distributions can be obtained with a higher degree of 
confidence than those in the relative age distributions, as exemplified by the smaller 
spread among the peak values of the mass distributions derived. This implies that mass 
determinations are mostly insensitive to the approach adopted. We reiterate that as 
extensive a wavelength coverage as possible is required to obtain robust and internally 
consistent age and mass estimates for the individual objects, with reasonable uncer- 
tainties. Finally, we conclude that the actual filter systems used for the observations 
should be used for constructing model colours, instead of using conversion equations, 
to achieve more accurate derivations of ages and masses. 

Key words: methods: data analysis - galaxies: spiral - galaxies: starburst - galaxies: 
star clusters 
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1 INTRODUCTION 

The systematic uncertainties in the determination of the 
ages, masses, and - to a lesser extent - extinction and metal- 
licity of young extragalactic star clusters, based on broad- 
band imaging observations, but using a variety of analy- 
sis approaches are poorly understood. Our aim in this pa- 
per is to intercompare the results obtained from Hubble 
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Space Telescope (HST) observations of carefully selected 
young star cluster samples using a variety of commonly used 
modelling techniques, and characterise the major differences 
among the techniques. We will do this specifically with re- 
spect to our own extensively tested multi-dimensional ap- 
proach, which we will use as our main benchmark. 

1.1 Star clusters as tracers of violent star forming 
episodes 

At one time or another during its lifetime, every galaxy will 
be affected by the external gravitational effects, however mi- 
nor, exerted by its neighbour galaxies. Irrespective of the 
precise interaction geometry, the results of such interactions 
are often striking: depending on the available gas reservoir, 
violent star formation will ensue, frequently predominantly 
in the guise of active star cluster formation (e.g., Whitmore 
ct al. 1999, de Grijs et al. 2001, 2003a,b, and references 
therein). Thus, where the above scenario holds, the age dis- 
tribution of a galaxy's star cluster population reflects its 
violent interaction history. 

The study of young and intermediate-age star cluster 
systems in a variety of galactic environments out to < 100 
Mpc has become a major field in extragalactic astrophysics 
in recent years, significantly stimulated by the superb imag- 
ing quality of the HST. One of the key diagnostic tools often 
utilised to predict the fate of such cluster systems is the clus- 
ter luminosity function (CLF; based on broad-band imaging 
observations) . 

Significant age spreads in young cluster systems - which 
might still be undergoing active cluster formation - affect 
the observed CLF (Meurer 1995, Fritze-v. Alvensleben 1998, 

1999, de Grijs et al. 2001, 2003b). This might, in fact, make 
an intrinsically log- normal CLF appear as a power- law (e.g., 
Miller et al. 1997, Fritze-v. Alvensleben 1998); the exact 
shape of the intrinsic CLF, whether a power law or a log- 
normal distribution, is still being debated (e.g., Vesperini 

2000, 2001, vs. Fall & Zhang 2001; see also Lamers et al. 
2004). It is obviously very important to obtain accurate age 
estimates for the individual clusters within such a system 
and to correct the observed CLF to a common age, before 
interpreting it as an intrinsic CLF (Fritze-v. Alvensleben 
1999, de Grijs et al. 2001, 2003b). 

1.2 Star cluster metallicities and the importance 
of interstellar dust 

The metallicities of star clusters produced in the high- 
pressure environments of galaxy interactions, mergers and 
starbursts are an important discriminator against the old 
Milky Way- type globular clusters (GCs), thought to be the 
oldest building blocks of most nearby spiral and elliptical 
galaxies. They are expected to correspond to the interstel- 
lar medium (ISM) abundances of the interacting/starburst 
galaxies, and are therefore most likely significantly more 
metal-rich than those of halo GCs in the Milky Way and 
other galaxies with old GC systems. However, ISM abun- 
dances span a considerable range among different galaxy 
types, from early- type spirals to dwarf irregulars (e.g., Fer- 
guson et al. 1998), and may also exhibit significant radial 
gradients (Oey & Kennicutt 1993, Zaritsky, Kennicutt & 



Huchra 1994, Richer & McCall 1995). Hence, a considerable 
metallicity range may be expected for star clusters produced 
in interactions of various types of galaxies and even among 
the clusters formed within one global galaxy-wide starburst. 

A significant increase of the ISM abundance in massive 
gas-rich galaxies may occur during strong bursts (Fritze- 
v. Alvensleben & Gerhardt 1994, their Fig. 12b). At the 
same time, typically within a few x 10 8 yr, some fraction of 
the gas enriched by dying first-generation burst stars may 
well be shock-compressed to cool fast enough to be built 
into later generations of stars or clusters produced in the 
ongoing burst. The same effect may occur when multiple 
bursts occur in a series of close encounters between two or 
more galaxies before their final merger. 

Precise (relative) metallicity determinations for individ- 
ual young star clusters are not only important to address 
these issues, but also for the correct derivation of ages from 
broad-band colours or spectral energy distributions (SEDs). 

Dust extinction is often very important in young cluster 
systems. In particular the youngest post-starburst galaxies 
and galaxies with ongoing starbursts often show strong and 
patchy dust structures. For instance, the youngest clusters 
in the overlap region of the two galactic discs in the Anten- 
nae galaxies are highly obscured in the optical and are best 
detected in near or mid-infrared observations (Mirabel et 
al. 1998, Mengel et al. 2001). Similarly, Lancon et al. (2003) 
discovered one of the most massive young star clusters in 
M82 based on near-infrared (NIR) HST observations; at op- 
tical wavelengths the cluster is completely obscured. Older 
merger remnants like NGC 7252 or NGC 3921 seem to have 
blown their inner regions clear of all the gas and dust left 
over from intense star formation (e.g., Schweizer et al. 1996). 
Extinction estimates toward individual clusters are therefore 
as important as individual metallicity estimates in order to 
obtain reliable ages and masses. 

1.3 Multi-passband photometry as a prime 
diagnostic 

Spectroscopy of individual massive young clusters, although 
feasible today with 8m-class telescopes for the nearest sys- 
tems, is very time-consuming, since observations of large 
numbers of clusters are required to obtain statistically sig- 
nificant results. Multi-passband imaging is a very interesting 
and useful alternative, in particular if it includes coverage 
of NIR and/or ultraviolet (UV) wavelengths (e.g., de Grijs 
et al. 2003c, Anders et al. 2004b). There are obviously lim- 
itations to the accuracy of the cluster parameters derived 
from broad-band imaging observations (e.g., de Grijs et al. 
2003b,c, Anders et al. 2004b, Bastian et al. 2004), but the 
relative overall characteristics derived for the cluster popu- 
lations as a whole appear to be relatively robust. 

In this paper we assess the systematic uncertainties in 
age and mass determinations, and to a lesser extent also 
in extinction and metallicity determinations, for young star 
cluster systems based on the use of broad-band, integrated 
colours, employing a variety of independently developed 
methods to analyse extragalactic star clusters as so-called 
"simple stellar populations" (SSPs): star clusters are the 
simplest objects to model, since they approximate single- 
age, single-metallicity populations with a range of stellar 
masses. Stochastic sampling effects of the stellar initial mass 
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function (IMF) also affect star cluster properties, in partic- 
ular for low-mass objects. However, since they affect broad- 
band photometry to a smaller extent than spectroscopy, and 
because we are dealing here with high-mass clusters only, we 
will not include these effects in this paper. Our main aim in 
this paper is to intercompare the results obtained for sets 
of well-calibrated cluster SEDs using a variety of commonly 
used modelling techniques, specifically with respect to our 
own extensively tested multi-dimensional approach (see Sec- 
tion 03 . 

In order to determine the absolute systematic uncer- 
tainties caused by the intrinsic differences in the models 
and methods in use in the literature, we distributed sets of 
broad-band star cluster photometry (described in Section^ 
to representatives of the various groups active in this field, 
and requested them to derive the cluster parameters using 
their specific methodology, wherever the data allowed this 
to be done. The models and methods are described in Sec- 
tion |3] we emphasize that most of the comparisons among 
methods done in this paper should be considered relative to 
the results obtained using the AnalySED method described 
in Section f3. 71 We summarise the results from applying our 
AnalySED approach to a large grid of artificial cluster data 
in Section |3] in order to establish the theoretical bench- 
mark for further comparisons among approaches. In Section 
01 we compare the overall, relative parameter distributions, 
while in Section|S]one-to-one comparisons between the vari- 
ous method+model combinations for the individual clusters 
in both of our samples are discussed; the results from each 
of the methods used for both cluster samples are included 
in Appendix A. We extend this discussion by considering 
the effects of converting the cluster photometry to different 
filter systems (Section|7J, and conclude the paper in Section 

El 



2 THE CLUSTER DATA SETS 

The field of stellar population modelling using extragalactic 
compact star cluster systems has undergone a major expan- 
sion since high-resolution, well-calibrated HST observations 
became available to the community. The application of stel- 
lar population synthesis of galactic subcomponents has be- 
come almost trivial for galaxies within ~ 20 — 30 Mpc, while 
star cluster population modelling is very well feasible out 
to ~ fOO Mpc, at least for the brighter (and therefore more 
massive) clusters within a given cluster population (see, e.g., 
de Grijs et al. [2003d] and Pasquali, de Grijs & Gallagher 
[2003] for examples toward and close to the distance limit). 
Ongoing and future HST programmes will continue to pro- 
vide high-resolution UV-optical-NIR imaging of large sam- 
ples of galaxies out to these distances. We therefore expect 
that modelling simple stellar populations and their broad- 
band SEDs will remain a key diagnostic tool for both the 
study of the evolution of extragalactic star cluster systems 
and their relation to Milky Way-type GCs, and for the anal- 
ysis of galactic star formation and interaction histories. 

Therefore, we decided to focus our comparison of mod- 
els and methods on HST-based imaging data. In addition, 
the calibration of HST measurements is well-understood 
and does, therefore, not introduce additional uncertainties 
as caused by, e.g., temporal variations in the atmospheric 



transmission that one would have to deal with if ground- 
based observations were used. 

We selected subsamples from large populations of young 
star cluster systems extensively studied in the literature, 
which we required to be among the brighter members of 
their respective populations (thus minimising the observa- 
tional uncertainties), as well as spanning a large age range 
(based on preliminary analyses, as described below). Ideally, 
we would have preferred to select cluster samples for which 
both HST measurements in a minimum of four broad-band 
passbands could be obtained, as well as independently de- 
termined parameters from spectroscopic observations. Un- 
fortunately, such data sets are yet not available, however. 
On the other hand, Schweizer, Seitzer & Brodie (2004) re- 
cently showed convincingly that spectroscopic age determi- 
nations are not necessarily better or more accurate than 
photometrically obtained ages, at least in the age range from 
~ 100 - 500 Myr. 

Our basic cluster samples were taken from the following 
sources: 

• NGC 3310, a nearby spiral galaxy exhibiting dominant 
star cluster formation in a circumnuclear starburst ring. This 
galaxy was covered by HST by the largest possible wavelenth 
range. The full set of eight broad-band STMAG magnitudes 
for the ~ 300 clusters in the galaxy's centre (located in the 
starburst ring and outside of it), from F300W ("mid-UV") 
to F205W were published and analysed in de Grijs et al. 
(2003a,c). We selected 17 of these clusters for the present 
analysis, all with well-determined, high-quality photometry 
in the entire set of available passbands (see Table 0. 

• NGC 4038/39 (the "Antennae"). Standard Johnson- 
Cousins UBVI photometry and Ha equivalent widths (EWs) 
for its large population of young to intermediate-age star 
clusters were obtained during a number of HST imaging 
campaigns by Whitmore and collaborators (see, e.g., Whit- 
more et al. 1999). Of the 20 objects selected for the pur- 
pose of this paper (Table 0, all have well-determined UBVI 
magnitudes (obtained via conversion of the HST flight sys- 
tem magnitudes using the Holtzman et al. [1995] conversion 
equations; but see Section 0, while 10 of them have mea- 
sured Ha EWs as well. 



3 MODELS AND METHODS 

Although the methods used to derive the global parameters 
of our cluster samples each have their own merits and dis- 
advantages, there is significant overlap among both the ex- 
tinction laws and the simple stellar population (SSP) models 
used for the stellar synthesis modelling. 

Therefore, we will first summarise the main character- 
istics of the SSP models and extinction laws used in this 
project. Subsequently, in Sections !3.3H3.7l each of the meth- 
ods employed to obtain the basic cluster parameters are de- 
scribed in detail, roughly in order of increasing complexity 
and sophistication. 

3.1 Extinction laws 

In Sections 13. 3H3.7I below, we will use a variety of Galactic 
extinction laws, as published by Savage & Mathis (1979; 
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Table 1. HST STMAG photometry of the NGC 3310 cluster sample 

The magnitudes are expressed in the STMAG HST flight system, derived from the count rates in the images as Tristmag — — 2.5 X 
log(counts s" 1 ) + zero-point offset; the zero-point offset is defined by the image header keywords PHOTFLAM and PHOTZPT. 



ID "1 F 300W "1F336W ™F439W »™F606W ™F814W ™F110W ™F160W "IF205W 
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17 
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18 


938 
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858 
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22 
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091 


22 


691 
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19 
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080 


19 


352 


± 





077 


19 


398 


± 





078 


20 


439 
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095 
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958 


± 





057 


21 


482 
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22 
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063 
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089 


± 





098 
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19 


500 


± 





095 


19 


858 
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744 
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762 
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22 


412 
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034 
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20 


032 
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174 
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002 
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740 
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079 
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311 
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069 
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377 
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079 


Gl 
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17 


523 
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023 


17 


043 
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18 


402 
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035 


18 


790 
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040 


20 


202 
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260 
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448 
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173 
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646 
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10 


788 
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011 


17 


248 
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030 
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411 
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007 
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870 
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124 
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668 
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134 
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493 
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25 


675 
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078 


19 


857 
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066 
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091 
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079 
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896 
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089 
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730 
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092 


22 


065 
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141 


22 


801 
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187 
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940 
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100 
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10 


18 


186 


± 





040 


18 


182 


± 





038 


18 


744 


± 





068 


19 


127 


± 





052 


20 


279 


± 





075 


20 


925 


± 





049 


21 


911 


± 





043 


22 


403 


± 





080 


Gl 


11 


19 


174 


± 





071 


19 


305 


± 





059 


19 


511 


± 





080 


20 


259 


± 





094 


20 


465 


± 





104 


21 


390 
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160 


21 


742 
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160 


22 


770 
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336 
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17 


407 
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17 


864 
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574 
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13 
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694 
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115 
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563 
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234 


19 


594 


± 
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057 
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199 
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284 
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029 
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098 
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955 
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072 


22 


847 


± 





107 
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14 
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022 
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370 
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219 
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219 
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308 
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250 
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032 
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175 
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208 
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109 
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774 
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105 


22 


101 
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094 


23 


289 
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217 
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870 
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19 


891 
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384 
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071 
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024 
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068 
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604 
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128 


± 





100 
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244 
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149 


± 
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063 
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762 
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22 


880 
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17 


19 


787 
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094 
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967 
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19 


455 
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080 
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818 
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117 


20 


903 
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080 


21 


218 
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056 


21 


009 


± 





056 


22 


234 


± 





059 



Table 2. Johnson-Cousins photometry and Ho EWs of the NGC 
4038/39 cluster sample 

The broad-band photometry is in magnitudes; lcr photometric 
uncertainties are on the order of 0.08 mag for all passbands. The 
Ha EWs are expressed in A. 



ID 




V 




B 




V 




J 


log( EW Hq ) 


G2-01 


21 


944 


21 


779 


21 


550 


21 


098 




G2-02 


21 


982 


21 


827 


21 


569 


21 


098 




G2-03 


23 


118 


23 


075 


22 


602 


21 


694 




G2-04 


20 


437 


21 


001 


20 


679 


20 


086 


2.942 


G2-05 


21 


786 


21 


870 


21 


212 


19 


995 


2.481 


G2-06 


21 


759 


21 


659 


21 


508 


21 


148 


0.700 


G2-07 


18 


407 


19 


145 


19 


010 


18 


634 


1.328 


G2-08 


21 


534 


21 


922 


21 


292 


20 


346 


2.853 


G2-09 


21 


388 


21 


412 


20 


800 


20 


041 




G2-10 


20 


273 


20 


732 


20 


290 


19 


516 


1.350 


G2-11 


23 


801 


23 


707 


22 


769 


21 


540 




G2-12 


18 


000 


18 


831 


18 


700 


18 


500 


2.487 


G2-13 


24 


430 


24 


271 


23 


316 


22 


181 




G2-14 


23 


377 


23 


073 


22 


201 


21 


157 


3.728 


G2-15 


18 


557 


19 


428 


19 


064 


18 


919 


3.497 


G2-16 


20 


159 


20 


511 


20 


321 


19 


823 


0.348 


G2-17 


19 


420 


19 


966 


19 


050 


18 


928 




G2-18 


22 


353 


22 


190 


21 


931 


21 


530 




G2-19 


23 


285 


23 


332 


22 


544 


21 


511 




G2-20 


24 


182 


23 


724 


22 


809 


21 


602 





Sections l3,4l and l3.6L Rieke & Lebofsky (1985; Section f3. 311 . 
Voshchinnikov & Il'in (1987; Section I3.4L and Fitzpatrick 
(1999; Section l3.5[l . as well as the starburst galaxy extinction 
law of Calzetti et al. (1994; Section Em. 

In the left-hand panel of Fig.^we show these extinction 
laws in relation to each other over the wavelength range of 
interest for the present study, normalised at an extinction of 
1 mag in the V band at 5500A, Av = 1 mag. In the right- 
hand panel we zoom in to display the differences among the 
individual extinction laws from 0.5 to 2.0 /im. From a com- 
parison of the individual extinction curves in the right-hand 
panel, it is clear that the differences are generally < 0.05 mag 
at wavelengths longward of 1/xm and shortward of ~ 0.8/im 
(with the exception of the Voshchinnikov & Il'in [1987] ex- 
tinction law). In the intermediate wavelength range, the dif- 
ferences are mainly driven by the Rieke & Lebofsky (1985) 
Galactic extinction law on the one hand and the Voshchin- 




A (/um) A (/Lim) 

Figure 1. Comparison of the various extinction curves adopted; 
the right-hand panel is a zoomed-in version of the left-hand 
panel. Galactic extinction laws: solid lines - Savage & Mathis 
(1979); dotted lines - Rieke & Lebofsky (1985); short-dashed lines 
- Voshchinnikov & Il'in (1987): long-dashed lines — Fitzpatrick 
(1999). Starburst galaxy extinction law: dot/long-dashed lines - 
Calzetti et al. (1994). 



nikov & Il'in (1987) curve on the other. Nevertheless, repre- 
sentative differences from the mean generally do not exceed 
0.1 mag, even at these wavelengths, and are often signifi- 
cantly smaller. 

If we place these differences among the extinction laws 
in their proper context, i.e., in comparison with the obser- 
vational data presented in Section |5J we see that they are 
generally of the same order or smaller than the observational 
uncertainties. In addition, as we will see in Section^ the vast 
majority of the sample clusters are characterised by Av <C 1 
mag, so that the differences among the various extinction 
laws become negligible for our sample clusters, irrespective 
of the analysis approach adopted. This conclusion is further 
strengthened if we realise that the wavelengths crucial for a 
successful determination of, in particular, extinction values 
and metallicities are the bluest optical/UV and the redder 
NIR passbands (see, e.g., de Grijs et al. 2003c, Anders et al. 
2004b), where the differences among the individual extinc- 
tion curves are smallest. 
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3.2 Simple stellar population models 

All of the methods used in this study rely on a comparison 
of the observational broad-band SEDs with a grid of model 
SEDs, in the sense that the star clusters are assumed to rep- 
resent "simple" stellar populations, i.e., single-age, single- 
metallicity populations with a range of stellar masses deter- 
mined by a given stellar IMF. The main differences among 
the various SSPs used for the comparison done in this paper 
are related to the use of different descriptions of the input 
physics, such as (i) a variety of stellar tracks or isochrones, 
which may or may not include critical phases in stellar evo- 
lution, such as the red supergiant (RSG) phase, Wolf-Rayet 
stars, and the thermally-pulsing AGB (TP-AGB) phase; (ii) 
slightly different IMF descriptions; and (iii) different (or 
no) treatment of nebular emission, particularly in the early 
phases of stellar evolution. 

Each of the approaches used in this paper is based on a 
comparison with a specific set of stellar evolutionary synthe- 
sis models. We will use (i) the most commonly used set of 
models developed by Bruzual & Chariot (1993, 1996 [BC96], 
2000 [BC00]; recently updated, 2003), (ii) the Starburst99 
models of Leitherer & Heckman (1995) and Leitherer et al. 
(1999), which are specifically matched to analyse the young 
stellar populations in starburst and interacting galaxies, and 
(iii) the Gottingen SSP models GALEV (Kurth et al. 1999, 
Schulz et al. 2002, Anders & Fritze-v. Alvensleben 2003). 



3.2.1 The Bruzual & Chariot SSP models 

The basic assumptions of the modern sets of the Bruzual 
& Chariot SSP models (used in Sections 13.31 13.51 and 13.61 . 
were first developed in Chariot & Bruzual (1991). The ver- 
sions of the code used in this paper, BC96 and BC00, in- 
clude a description of SSP evolution for a range of metallic- 
ities from Z = 0.0004 to Z = 0.10. The BC96 models and 
more recent versions are (mostly) based on the evolutionary 
tracks of the Padova group (Bressan et al. 1993, Fagotto et 
al. 1994a, b,c; with additional empirical spectra for stellar 
masses between 0.1 and 0.7 Mq; Chariot & Bruzual 1991), 
and cover an age range from 1.25 x 10 5 to 2 x 10 10 yr, typ- 
ically computed using 220 unequally spaced age intervals. 
The important TP-AGB phase (see Section r3.2.3H is treated 
semi-empirically. These models cover all of the important 
phases of stellar evolution from the zero-age main sequence 
to the post- AGB phase and beyond, for stars with effective 
temperatures of 2000 < T cff (K) < 50,000. The stellar spec- 
tra are based on the theoretical spectral library compiled by 
Lejeune, Cuisinier & Buser (1997, 1998), which use the the- 
oretical stellar atmosphere calculations of Kurucz (1979), 
Fluks et al. (1994), and Bessell et al. (1989, 1991). Leje- 
une et al. (1997, 1998) have corrected the stellar continuum 
shapes from the Kurucz (1979) models to agree with ob- 
served colours from the UV to the K band. SSPs, covering 
the wavelength range from the extreme UV (5A) to the far- 
infrared (100/im) with a resolution depending on the spec- 
tral range, were calculated assuming a Salpeter (1955)-type 
IMF, £(m) oc m~ a with a = 2.35 and masses ranging from 
~ 0.1M Q up to 125 M . 



3.2.2 The Starburst99 models 

The Leitherer et al. (1999) Starburst99 models (used in Sec- 
tions 13.41 and l3,bl constitute an improved and extended ver- 
sion of the suite of models initially published by Leitherer 
& Heckman (1995). These models were specifically devel- 
oped for the evolutionary synthesis analysis of populations 
of massive stars, and are best suited to the conditions typ- 
ically found in starburst environments. They are based on 
the Geneva stellar evolution models and the new model at- 
mosphere grid compiled by Lejeune et al. (1997). The tracks 
of Meynet et al. (1994) were used for stars with masses in 
excess of 12-25 Mq (depending on metallicity), with the en- 
hanced mass loss prescription in order to better approximate 
most Wolf-Rayet properties (except the Wolf-Rayet mass- 
loss rate itself) compared to the standard mass loss scenario. 
For stars with masses in the range 0.8 < m*/M0 < 12, they 
used the standard mass- loss tracks of Schaller et al. (1992), 
Schaerer et al. (1993a,b) and Charbonnel et al. (1993). These 
tracks include the early AGB evolution until the first ther- 
mal pulse for stars with masses m* > 1.7M0. The Star- 
burst99 models also include observational high-resolution 
UV spectra, to allow for the analysis of stellar and interstel- 
lar absorption lines and line profiles at various metallicities. 

The SSP models cover an age range between 10 6 and 
10 9 yr, with an age resolution of 0.1 Myr, for all five metal- 
licities, Z = 0.001, 0.004, 0.008, 0.02, and 0.04 over the entire 
spectral range from the extreme UV to the infrared. Nebular 
continuum emission is included in the models in a simpli- 
fied fashion; its contribution becomes important when hot 
stars providing ionising photons (and thus line emission) are 
present (see Section 13.2.31 . 

The synthesised models that we use in Section ^. 6l below 
were calculated using a standard Salpeter-type IMF, char- 
acterised by stellar masses in the range 1 < m^/M© < 100. 

3.2.3 The GALEV Gottingen SSP models 

The GALEV SSPs (used in Sections 0O1 EH and 03 are 

based on the set of stellar evolutionary tracks (Kurth et 
al. 1999), and in later versions the isochrones (Schulz et al. 
2002, Anders & Fritze-v. Alvensleben 2003), of the Padova 
group (with the most recent versions using the updated 
Bertelli et al. [1994; and unpublished] isochrones; the latter 
also include the TP-AGB phase) for the metallicity range of 
0.0001 < Z < 0.05, tabulated as five discrete metallicities 
(Z = 0.0004, 0.004, 0.008, 0.02, and 0.05, corresponding to 
[M/H] « [Fe/H] = -1.7,-0.7,-0.4,0.0, and +0.4, respec- 
tively). For lower-mass stars (0.08 < m,/Mg < 0.5), which 
contribute very little to the integrated light of young and 
intermediate-age SSPs governed by any standard, Salpeter- 
type IMF, the Padova models are supplemented with the 
Chabrier & Baraffe (1997) theoretical calculations that in- 
clude a new description of stellar interiors of low-mass ob- 
jects and use non-grey atmosphere models. 

The GALEV models are furthermore once again based 
on the theoretical stellar libraries of Lejeune et al. (1997, 
1998) for a broad range of metallicities. For stars hotter 
than T e ff = 50, 000 K, pure black-body spectra are adopted, 
as for the Bruzual & Chariot models. The full set of models 
spans the wavelength range from 90A to 160/im. 

The Salpeter-type IMF assumed is characterised by a 
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lower cut-off mass of 0.15Mq; the upper-mass cut-off ranges 
between 50 and 70 M Q , and is determined by the mass cov- 
erage of the Padova isochrones for a given metallicity. 

Kurth ct al. (1999) cover ages between 1 x 10 7 to 
1.6 x 10 10 yr, with an age resolution of 10 7 , 10 8 , and 10 9 
yr, for ages < 10 s , between 10 s and 10 9 and > 10 9 yr, re- 
spectively. Schulz et al. (2002) and Anders & Fritze-v. Al- 
vensleben (2003) extended the age range down to 4 x 10 6 yr 
(and slightly reduced the upper age limit to 14 Gyr), while 
improving the age resolution to 4 Myr for ages up to 2.35 
Gyr, and 20 Myr for greater ages. 

The Schulz et al. (2002) version includes important im- 
provements with respect to the older versions; they use the 
newer Padova isochrones which include the important stellar 
evolutionary TP-AGB phase. At ages ranging from ~ 100 
Myr to - 1 Gyr, TP-AGB stars account for 25 to 40 per 
cent of the bolometric light, and for 50 to 60 per cent of the 
K-b&nd emission of SSPs (see Chariot 1996, Schulz et al. 
2002). Schulz et al. (2002) show that the effect of including 
the TP-AGB phase results in redder colours for SSPs with 
ages between ~ 10 8 and 10 9 yr, with the strongest effect (up 
to > 1 mag) being seen in (V — K) for solar metallicity, and 
in [V — I) for Z > 0.5Zq. Shorter-wavelength colours and 
lower metallicity SSPs are less affected. Since most young 
to intermediate-age star cluster systems observed in HST 
passbands equivalent to the standard V and I filters, are in 
fact aged between about 100 Myr and 1 Gyr, and have often 
close-to solar metallicities, inclusion of the TP-AGB phase 
in the models is obviously important. 

Finally, Anders & Fritze-v. Alvensleben (2003) included 
gaseous continuum emission and an exhaustive set of nebu- 
lar emission lines to the GALEV suite, assuming comparable 
metallicities for the star cluster and the surrounding ionised 
gas. Nebular emission is shown to be an important contribu- 
tor to broad-band fluxes during the first few x 10 7 yr of SSP 
evolution, the exact details depending on the metallicity. 



3.2.4 Model comparison 

First, we present a basic comparison among the SSP mod- 
els used in this paper. Kurth et al. (1999) and Schulz et al. 
(2002) concluded that, compared to the models of BC96 and 
Bruzual & Chariot (1993), respectively, their sets of galev 
models agree very well for solar metallicity and a Salpeter 
IMF, for (B — V) colours, and from the UV up to ~ 7000 A, 
respectively. However, between 7000A and 12,000A, as well 
as in the NIR H and K band regime, the BC96 flux con- 
tribution is considerably lower than that of the Schulz et 
al. (2002) spectrum at those wavelengths, which they at- 
tribute to the different treatment of the TP-AGB evolution- 
ary phase (with T aS ~ 3000 K). 

Anders & Fritze-v. Alvensleben (2003) compare the 
most up-to-date galev models that include nebular line and 
continuum emission with the Starburst99 models. Despite 
the differences in the input physics and the different sets of 
stellar tracks and/or isochrones used by these teams, they 
conclude that the differences between both sets of models are 
minor at short optical wavelengths (e.g., (B — V) colours) 
during the first Gyr of evolution, and are mainly due to the 
better time resolution of the Starburst99 models. Longer- 
wavelength comparisons show larger differences, due to the 



different input physics, and in particular a different treat- 
ment of RSGs and the TP-AGB phase. 

In summary, it appears that over most of the optical 
wavelength range all of the commonly used SSP models are 
fairly similar, with minor differences depending on the de- 
tailed input physics and the treatment of the various evolu- 
tionary phases. At longer (NIR) and shorter (bluer) wave- 
lengths, the differences become more significant, and will 
lead to systematic differences in the determination of the 
basic properties of SSPs, as we will see below. 

Quantitatively, for a given set of input physics, varying 
parameters including the IMF slope, mass loss and convec- 
tion prescriptions, one can justify differences of up to ~ 0.05 
mag in (B — V) (e.g., Yi 2003). However, the difference be- 
tween - for instance - the Padova and Geneva stellar evo- 
lutionary tracks when used by an otherwise identical SSP 
code amounts to A(B — V) ~ 0.05 mag, and even more in 
A(V - I) and A(V - K) (e.g., Leitherer et al. 1996, Schulz 
et al. 2002). 

It is clear from the outset that the application of the 
various SSP models to our sets of sample clusters will re- 
sult in significantly different masses and mass distributions, 
simply because of the different low and high-mass bound- 
aries adopted for the Salpeter-type IMF. The mass ratios 
expected to result from the Starburst99 : the Bruzual & 
Chariot : the GALEV SSPs are 1 : 22.4 : 33.1, or in loga- 
rithmic mass units, mass estimates based on the Bruzual & 
Chariot SSPs will result in masses that are 1.35 dex higher 
than those from the Starburst99 models; the GALEV masses 
will be 1.52 dex more massive than the Starburst99 ones. 



3.3 Method 1: Optical/NIR sequential analysis 
("Sequential O/IR") 

The Sequential O/IR method is a two-step approach to 
derive the age, metallicity and extinction values associ- 
ated with a given cluster. First, the extinction is estimated 
using the BVI passband combination; subsequently, the 
extinction-corrected, intrinsic colours are compared, in a 
least-squares sense, to the BC96 SSP models in order to 
estimate the cluster age. 

While the (B — V) vs. (V — I) colour-colour diagram is 
affected by the well-known age-metallicity degeneracy, the 
age and extinction trajectories are not entirely degenerate 
for this particular choice of optical colours. For SSPs older 
than ~ 100 Myr (i.e., (B - V) > 0.18 mag), all age trajec- 
tories show the same, roughly linear growth of the (V — I)o 
vs (B — V)o colours, irrespective of their metallicity. As a 
consequence, for such ages, (B — V) vs. (V — I) SSP analy- 
sis enables us to derive the visual extinction (and therefore 
the intrinsic colours and magnitudes), prior to any age and 
metallicity estimates (cf. de Grijs et al. 2001). 

Using the intrinsic colours we can now derive the most 
representative age (and metallicity), by minimising (in a 
least-squares sense), 



CF mtr -Ch SSP (t,Z) 



(1) 



where Cl| ntr and Clf SP (t, Z) are the intrinsic and the model- 
predicted colour indices in a given colour denoted by i, re- 
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Figure 2. (V — I)o vs. (V — J)o colour-colour diagram of the 
BC96 SSPs. The thin curves represent isometallicity tracks; the 
thick curves connecting filled symbols represent isochrones (from 
left to right, log t/yr = 8.5, 9.0, 9.5 and 10.0, respectively). As 
an example, the thin boxes represent the medians of typical ltr 
error distributions for SSPs of solar metallicity at ages of 600 Myr 
and 10 Gyr. The arrow indicates the effects of reddening on these 
models, for a visual extinction Ay = 1 mag. For the sake of clarity, 
only SSPs older than 100 Myr are shown. 



spectively, for SSPs with age t and metallicity Z; aci; are 
the lcr uncertainties. 

If NIR photometry is available, the cluster age and 
metallicity can be derived simultaneously: for instance, the 
isochrones and iso-metallicity tracks define a grid in (V — I) 
vs. (V — J) space, thus allowing one to lift the age-metallicity 
degeneracy, as shown in Fig. [5] Finally, cluster masses are 
obtained from their luminosities via the age and metallicity 
dependent mass-to- light ratio. 

Our recent study of the intermediate-age star cluster 
population in M82's region B provides a good example of 
what can be achieved if both optical and NIR data are 
available. In that case, both the age-extinction and the age- 
metallicity degeneracies can be lifted. Further details, and 
a discussion about how photometric errors propagate into 
extinction, age, metallicity and mass uncertainties are given 
in Parmentier, de Grijs & Gilmore (2003). 



3.4 Method 2: Reddening-free Q parameter 
analysis ("Q-Q") 

The basic Q-parameter analysis is a powerful method to de- 
termine SSP ages and extinction values independently. The 
host galaxy's internal extinction to a given cluster is derived 
from the Qubv parameter (Johnson & Morgan 1953), 



Qubv = (U 
= (U- 



B) -[E(U-B)/E(B- 
B) — 0.72 x (B — V), 



V)} x(B- V) 
(2) 



To estimate the internal reddening we assess the loci of the 
clusters in the (B — V) vs Qubv plane, compared to the 
intrinsic (B — V)o colours of the Starburst99 and GALEV 
SSPs. 

Subsequently, age estimates are obtained by minimis- 
ing the clusters' loci in the plane of the reddening-free Qubv 
and Qubv i parameters, where the latter is defined as (Whit- 



more et al. 1999) 
Qubvi = (U - 

= (u- 



B) - [E(U-B)/E(V 
B) — 0.58 x (V — I), 



I)] x (V - J) 
(3) 



with respect to the Starburst99 models. Unfortunately, due 
to a loop of the evolutionary tracks in Q-Q space, one cannot 
achieve accurate age estimates in the range of log(Age/yr) 
from ~ 6.5 to 7.2 (see Whitmore et al. 1999). The availability 
of Ha observations will greatly facilitate our age estimates 
in this age range. 

3.5 Method 3: Ha luminosities in addition to 
broad-band fluxes ("BB+Ha") 

For those clusters for which we have Ha flux or EW infor- 
mation available, we compared the five observed magnitudes 
(UBVIRa) with the predicted SEDs from the BC00 SSPs 
at solar metallicity. For each available model age, we var- 
ied the reddening between E(B — V) = 0.0 and 3.0 in steps 
of 0.02, and assume the Galactic extinction curve of Fitz- 
patrick (1999). Each model-age/reddening combination is 
scaled to match the cluster V-band magnitude, and then 
compared with the observations using a standard \ 2 min- 
imisation technique: 



2 

Xmin 



t,E(B - V)) = mini J2 



/ model 



obs\ 2 

m x ) 



(4) 



where A = U,B,V,I,Ha and w x = [(0.05) 2 + ir^]" 1 . Here, 
(Jm x is the photometric uncertainty (in magnitudes) for a 
given bandpass; we have included an additional uncertainty 
of 0.05 mag, which represents the uncertainties in the models 
themselves (Section 13.2.41 see Fall, Chandar & Whitmore 
[2004] for validation of the method). Since bright objects 
can have very small photometric uncertainties, this keeps 
the weights from "blowing up" in the very bright object 
regime. 

The predicted Ha model flux is calculated from the 
total number of ionising photons under the assumption of 
photon-limited case B recombination. When converting to 
magnitudes, we determined the zero-point offset between the 
models and observations empirically, by comparing the ob- 
served difference in (mu a — my) for the strongest Ha emit- 
ters with model predictions. We then applied the empirically 
determined zero-point offset to the entire dataset, and found 
that for clusters younger than 10 Myr, with measurable Ha 
emission, age estimates were in good agreement with those 
derived by Whitmore & Zhang (2002). (The measured Ha 
fluxes for the clusters were converted to the VEGAMAG sys- 
tem using the prescription given in the WFPC2 Data Hand- 
book.) 

We note that for those clusters without measurable Ha 
EWs, we applied the simplified method based on the broad- 
band luminosities only, but using otherwise the same proce- 
dure. 



3.6 Method 4: Three-dimensional SED analysis 
("3DEF") 

The next step up in complexity of the fitting algorithms 
used involves the fitting of the observed cluster SEDs to 
the Starburst99 and BC00 models using a three-dimensional 
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(3D) maximum likelihood method, 3DEF, with respect to a 
pre-computed grid of SSP models. This procedure was de- 
scribed in detail by Bik et al. (2003) , based on their analysis 
of archival HST observations of the central star clusters in 
M51, and applied successfully to the intermediate-age star 
cluster system in M82 B by de Grijs et al. (2003b) and to the 
extended cluster sample in M51 by Bastian et al. (2004). The 
initial cluster mass Mi, age and extinction E(B — V) were 
adopted as free parameters. For those clusters with upper 
limits in one or more filters, but still leaving us with a min- 
imum of three reliable photometric measurements, we use 
a two-dimensional maximum likelihood fit ( "2DEF" ) , using 
the extinction probability distribution for E(_B — V). This 
distribution was derived for the clusters with well-defined 
SEDs over the full wavelength range (see Bik et al. 2003 for 
a full overview of this procedure) . The derivation of the most 
representative set of models for a given cluster is done via 
a least-squares (x 2 ) minimisation technique, in which the 
observed cluster SED is compared to the full grid of SSP 
models. In the application of the 3DEF-method Bik et al. 
(2003) and Bastian et al. (2004) assumed an uncertainty of 
0.05 mag in the magnitudes of the cluster models (0.1 mag 
in the UV filters). 

3.7 Method 5: Multidimensional SED analysis 
("AnalySED") 

Finally, we have developed a sophisticated SED analysis tool 
that can be applied to photometric measurements in a given 
number N(N > 4) of broad-band passbands (see de Grijs 
et al. 2003c, Anders et al. 2004b). We apply a 3D \ 2 min- 
imisation to the SEDs of our star clusters with respect to 
the GALEV SSP models, to obtain the most likely combina- 
tion of age t, metallicity Z and internal extinction E(B — V) 
for each object (see Anders et al. 2004b; Galactic foreground 
extinction is taken from Schlegel, Finkbeiner & Davis 1998) . 

In order to obtain useful results for all of our three free 
parameters, i.e., age, metallicity and extinction 1 , we need a 
minimum SED coverage of four passbands. 

Each of the models is assigned a probability, determined 
by a likelihood estimator of the form p ~ exp(— x 2 ), where 

2 t < ryjj/r, xr\ \ \ " (m bs — ?ri mo dcl) 2 

X [t, Z,E(B - V),m c i) = - j '—. (5) 

models obs 

Clusters with unusually large "best" \ 2 are rejected, 
since this is an indication of calibration errors, features not 
included in the models (such as Wolf-Rayet star dominated 
spectra, objects younger than 4 Myr, etc.) or problems due 
to the limited parameter resolutions. We include an addi- 
tional 0.1 mag per passband for "model uncertainties" (0.2 
mag for UV filters; see Section [3.2.411 . 

Subsequently, the model with the highest probability is 
chosen as the "best-fitting model". Models with decreasing 
probabilities are summed up until reaching 68.26 per cent to- 
tal probability (i.e., the 1 a confidence interval) to estimate 

1 Strictly speaking, the cluster mass is also a free parameter. 
Our model SEDs are calculated for SSPs with initial masses of 
1.6x 10 9 Mq; to obtain the actual cluster mass, we scale the model 
SED to match the observed cluster SED using a single scale factor. 
This scale factor is then converted into a cluster mass. 



the uncertainties on the best-fitting model parameters. For 
each of these best-fitting models the product of the relative 
uncertainties (^S£_ x ma33+ x fi) was calculated (the su- 

v age mass Z ' * 

perscripts indicate the upper ( + ) and the lower limits (~), 
respectively). The relative uncertainty of the extinction was 
not taken into account, since the lower extinction limit is of- 
ten zero. For each cluster, the data set with the lowest value 
of this product was adopted as the most representative set of 
parameters. In cases where the analysis converged to a sin- 
gle model, a generic uncertainty of 30 per cent was assumed 
for all parameters in linear space, corresponding to an un- 
certainty of to' 15 dex in logarithmic parameter space. See 
also de Grijs et al. (2003a, b,c) and Anders et al. (2004a) for 
applications of this algorithm to NGC 3310 and NGC 6745, 
and NGC 1569, respectively, and Anders et al. (2004b) for 
a theoretical analysis of its reliability. 

We caution that the multi-passband combinations must 
not be biased to contain mainly short wavelength, nor 
mainly long- wavelength filters. Coverage of the entire op- 
tical wavelength range, if possible with the addition of UV 
and NIR data, is most preferable (de Grijs et al. 2003c, An- 
ders et al. 2004b). 

Finally, we emphasize once again that we will use this 
AnalySED method as the basis for our comparisons among 
the different approaches employed in this paper. This deci- 
sion is based on the fact that the method was validated and 
tested extensively, both empirically (de Grijs et al. 2003a,c) 
and theoretically (Anders et al. 2004b), so that we under- 
stand the systematic uncertainties inherent to this approach 
in depth. In the following section, we will summarise the re- 
sults from our extensive validation of the AnalySED method, 
in order to justify its use as our benchmark approach for 
comparison with the other methods descirbed in the previ- 
ous sections in the remainder of this paper. 



4 ESTABLISHING OUR BENCHMARK 
APPROACH WITH ARTIFICIAL DATA 

In Anders et al. (2004b) we presented a detailed study of the 
reliability and limitations of our AnalySED approach. We 
computed a large grid of broad-band HST-based star cluster 
SEDs on the basis of our GALEV models for SSPs, including 
all relevant up-to-date input physics for stellar ages > 4 Myr. 
We constructed numerous artificial cluster SEDs, and varied 
each of the input parameters (specifically, age, metallicity, 
and internal extinction; see Section 13.71 in turn to assess 
their effects on the robustness of our parameter recovery. 
For each clean model artificial cluster SED we calculated 
10,000 additional clusters, with errors distributed around 
the input magnitudes in a Gaussian fashion. 

By analysing artificial clusters, using a variety of input 
parameters with our AnalySED approach, we found in gen- 
eral good agreement between the recovered and the input 
parameters for ages < 10 9 yr, i.e., exactly the age range 
of interest for the clusters in NGC 3310 and the Antennae 
galaxies analysed in this paper. 

We considered several a priori restrictions of the full 
parameter space, both to the (correct) input values and to 
some commonly assumed values. We easily recover all re- 
maining input values correctly if one of them is restricted 



Systematic uncertainties in cluster parameters based on broad-band imaging 9 



T*e — 1 1 j. 

1 2 3 

cluster number 



□■o 



Z free - 
Z=0.008 *- 



Z=0.02=solar o 
input values 



10.5 
10 

J- 9.5 
o> 

S. 9 



5> 



8.5 



7.5 



6.5 



[] |T 

1 2 3 

cluster number 



Z free - 
Z=0.008 - 



Z-0.02-solar o 
input values 



7 
6.8 
6.6 
6.4 
6.2 



5.8 
5.6 



5.4 



111® 



2 3 
cluster number 




2 3 
cluster number 



Figure 3. Accuracy of the retrieval of the ages and masses of 
artificial star clusters based on a wavelength coverage as avail- 
able for our NGC 3310 cluster sample. The artifical clusters are 
characterised by E(B — V) = 0.1 mag, Z = Zq, and ages of 8, 
60, and 200 Myr, and 1 Gyr (objects 1-4). The different sym- 
bols represent the retrieved values based on a variety of a priori 
assumptions on the clusters' metallicities. 



a priori to its correct input value; this also provides a san- 
ity check for the reliability of our code. We conclude that 
the age-metallicity degeneracy is responsible for some mis- 
interpretations of clusters younger than ~ 200 Myr. If we 
restrict one or more of our input parameters a priori to in- 
correct values (such as by using, e.g., only solar metallicity, 
as often done in the literature), large uncertainties result in 
the remaining parameters. 

In order to provide a robust theoretical benchmark 
for the observational study of systematic uncertainties pre- 
sented in the remainder of this paper, here we re-validate 
our AnalySED approach using the exact filter combinations 
available for our NGC 3310 and Antennae cluster samples 
(Tables Q and , i.e., by computing a large grid of broad- 
band star cluster SEDs in a similar fashion as done in An- 
ders et al. (2004b) - although using only 1,000 additional 
artificial clusters to quantify our model uncertainties; the 
difference between this approach and the 10,000 additional 
artificial clusters used in Anders et al. (2004b) is negligi- 



Figure 4. Accuracy of the retrieval of the ages and masses of 
artifical star clusters based on a wavelength coverage as available 
for our Antennae cluster sample; technical details are as in Fig. 

m 



ble, however. Our artifical cluster SEDs were computed for 
ages of 8, 60, and 200 Myr and 1 Gyr (the age range cov- 
ered by our sample clusters; see Sections 15.11 and 15 . 21 . and 
for a fixed (internal) extinction of E(B — V) = 0.10 mag 
and solar metallicity. The extinction value adopted roughly 
corresponds to the mean extinction derived for the indi- 
vidual clusters (Sections 15.11 and 15 . 2fl : extinction variations 
among the sample clusters are small, and their effects (for 
the derived range of extinction values) on the age and mass 
estimates are negligible (de Grijs et al. 2003a, Anders et 
al. 2004b). The adopted solar metallicity also corresponds 
roughly to the mean metallicity derived for the sample clus- 
ters, although it may not be correct in individual cases. In 
order to quantify the effects of metallicity variations, we 
attempted to retrieve the ages and masses of our artificial 
clusters by assuming both the correct (solar) and incorrect 
(Z — 0.008 = QAZq ) metallicities as a priori restrictions. In 
addition, we retrieved the ages and masses of the artificial 
clusters without any restriction to the resulting metallicity 
("Z free"). The latter provides a quantitative indication of 
the importance of the age-metallicity degeneracy. 
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The results of this re-validation are shown in Figs. [3] 
and 2] for filter coverage as for NGC 3310 and the Anten- 
nae galaxies, respectively. It is immediately clear that the 
accuracy of the parameter retrieval is significantly better 
for the NGC 3310 clusters than for those in the Anten- 
nae galaxies, which simply reflects the available filter sets 
(cf. de Grijs et al. 2003a,c, Anders et al. 2004b). Never- 
theless, in all cases the ages are retrieved well within the 
modelling uncertainties, for any assumption on the clus- 
ters metallicity. For the coverage corresponding to the NGC 
3310 clusters, the difference between input and retrieved 
ages is Alog(Age/yr) < 0.3, and in the majority of cases 
Alog(Age/yr) < 0.15. Except for the 60 Myr-old cluster, 
Alog(Age/yr) < 0.3 also for those clusters covered by the 
same passbands as the Antennae clusters (although the age- 
metallicity degeneracy is somewhat more important for the 1 
Gyr-old cluster in this case). The corresponding uncertainty 
in the retrieved age of the 60 Myr-old artificial cluster is 
about twice as large as for the other clusters; its cause is 
unclear, since the retrieved extinction and metallicity val- 
ues for this object are not significantly more uncertain than 
for the other objects. 

A similar behaviour, i.e., with slightly larger uncertain- 
ties for the Antennae-equivalent wavelength coverage com- 
pared to the coverage of the NGC 3310 sample, is seen for the 
retrieved masses of the artificial clusters, although to a lesser 
extent. The mass uncertainties for all clusters older than 60 
Myr are Alog(M c i/M ) < 0.10 - 0.15. For the youngest, 8 
Myr-old object, we are only able to retrieve the masses to 
within several 0.1 dex in mass (somewhat more accurately 
for the NGC 3310-equivalent wavelength coverage, particu- 
larly if the adopted metallicity is close to the actual value); 
this is most likely caused by the uncertainties inherent to 
our present knowledge of stellar evolution in this age range 
(such as, e.g., the importance of the RSG phase), and the 
relatively coarse age resolution compared to the rapidity of 
changes in stellar evolution around 6-12 Myr. 

Thus we have shown, based on well-understood artificial 
data, that we understand quantitatively the uncertainties in- 
herent to using our AnalySED approach for age and mass 
determinations of star clusters based on broad-band imag- 
ing. We will be using this approach as our benchmark for 
comparing our results to those obtained using alternative 
methods commonly in use in the community. 

In Section 13.2.41 we concluded that the differences 
among the various prescriptions used for the input physics 
in modern sets of SSP models are very small indeed and ap- 
parently not biased systematically. As a consequence of the 
analysis performed in this section, we conclude then that 
any differences in the individual (as well as in the mean) 
cluster ages and masses that we will find in subsequent sec- 
tions (over and above the modelling uncertainties quantified 
here) are most likely caused by intrinsic differences among 
the various methods. 



5 COMPARISON OF THE RELATIVE AGE 
AND MASS DISTRIBUTIONS 

5.1 Extensive wavelength coverage: NGC 3310 

To start our comparison of methods, we will focus on the 
extensive wavelength coverage of the NGC 3310 star cluster 



system. With coverage from the F300W HST mid-UV pass- 
band to the NIR F205W passband, the resulting broad-band 
SEDs were shown to have sufficient leverage to distinguish 
metallicity, extinction, and stellar population (age) effects 
(de Grijs et al. 2003c). 

While a wavelength coverage as extensive as possible is 
preferred, the use of iTST-flight system magnitudes (cf. Sec- 
tion|^J limits the application of the NGC 3310 comparison to 
the use of the GALEV SSP models (see Section f3. 2. 31 . which 
we folded through the HST/WFPC2 filter curves ourselves 
(e.g., de Grijs et al. 2003c, Bastian et al. 2004). 

We emphasize, however, that we prefer to use the 
original HST-flight filter system, rather than conversions 
to "standard" systems; in Section [7] we will discuss the 
systematic effects unavoidably introduced when converting 
HST-flight system magnitudes to the "standard" Johnson- 
Cousins system. 

Figure|S|shows the resulting age distributions for the 17 
NGC 3310 star clusters used for this exercise, obtained using 
a variety of approaches. In Fig. [SJi, we display the relative 
age distribution of the NGC 3310 clusters based on the full 
multi-dimensional SED analysis (Section 13.71 . in which we 
left all of the cluster ages, masses, metallicities and extinc- 
tion values as free parameters. The vertical dotted line at 
log(Age/yr) = 6.6 denotes the lower age limit of the GALEV 
SSP models; the vertical error bars indicate the Poissonian 
uncertainties. 

Panel (b) provides a direct comparison of the effects of 
metallicity variations. Here, as well as in the other panels 
in this figure, we have adopted solar metallicity for the in- 
dividual star clusters. While this is not necessarily correct 
in general, restricting the metallicity to the solar value al- 
lows a more robust comparison among the various models 
and methods 2 . The differences in the relative age distribu- 
tions between Figs. Hjji and b are therefore entirely and ex- 
clusively due to the different assumptions on the clusters' 
metallicities. They reflect the well-known effects of the age- 
metallicity degeneracy (e.g., Ferreras & Yi 2004; see de Grijs 
et al. [2003a,c] for detailed studies of this effect in NGC 3310 
and NGC 6745). Based on the multi-dimensional analysis 
presented in Fig. HJ,, most (70 per cent) of the NGC 3310 
clusters in our current sample are characterised by metallici- 
ties Z < 0.01, with the remainder split evenly between solar 
metallicity (Z = 0.02) and Z ~ 0.05, see Fig. 03 Support for 
these metallicity estimates is provided by the unusually low 
(subsolar) metallicity found independently in star-forming 
regions surrounding the nucleus of NGC 3310, while the nu- 
cleus itself appears to have solar metallicity (e.g., Heckman 
& Balick 1980, Puxley, Hawarden & Mountain 1990, Pas- 
toriza et al. 1993). Similarly, the age-extinction degeneracy 
may contribute to some extent, although its effect is likely 
less than that of the age-metallicity degeneracy (cf. de Grijs 
et al. 2003c). The extinction in this sample of NGC 3310 

2 This is because by treating metallicity as a fit parameter we 
introduce the complexity of an additional free parameter, which 
has the potential to render the computational solution less stable 
and robust. In view of the small effects associated with small 
differences in metallicity, for the purpose of this exercise, we opt 
for the more robust approach to adopt a single (solar) metallicity. 
We note that this reflects common practice in the literature (but 
see de Grijs et al. 2003c). 
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Figure 5. Metallicity estimates of the NGC 3310 and Anten- 
nae sample clusters, based on the results obtained from the 
AnalySED multi-dimensional approach. Note that these arc first- 
order estimates (yielding an internally consistent overall metal- 
licity distribution; see Anders et al. 2004b), and that the un- 
certainties in the individual cluster metallicities are on the or- 
der of one step in our metallicity grid, covering metallicities of 
O.OO4Z Q ,O.2Z ,O.4Z Q ,Z (Z = 0.02), and 2.5Z e . 



Table 3. Characteristics of the overall relative mass distributions 
of our NGC 3310 cluster sample 
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clusters, E(£> — V), decreases with age - although with a 
large scatter - from E(B - V) ~ 0.35-0.40 mag at ~ 10 6 ' 5 
yr, to E(B - V) < 0.1 mag at ~ 10 8 yr. 

Figure [(J can be compared directly with Fig. |SJd; the 
only difference between these two panels is that we used 
the "3DEF" method (SectionOHD 1 instead of the AnalySED 
multi-dimensional approach. We used the GALEV SSP mod- 
els in both cases. It is encouraging to see that the use of 
either method results in very similar relative age distribu- 
tions. More quantitatively, the two peaks in the age distri- 
bution are reproduced to a very high degree of confidence 
(> 99 per cent), although their relative amplitudes are sub- 
ject to small-number statistics; as a result, the straightfor- 
ward application of a Kolmogorov-Smirnov (KS) test 3 yields 
a probability that these data points were drawn from sta- 
tistically different distributions of ~ 80 per cent, but with a 
very large uncertainty because of the small number of data 
points used. 

Finally, in panel (d) we have replaced the GALEV SSP 
models by the Starburst99 models (Section 13.2.20 . which re- 
sults in a markedly different age distribution. This is most 
likely caused by two effects, which are in essence the most 

3 From a pure statistics approach, the application of a KS test to 
our results is strictly speaking invalid. While the relatively small 
number of data points is in principle acceptable and sufficient 
for our purpose, systematic effects related to the sample selection 
are not taken into account, while they may in fact dominate. This 
thus renders the results not very illustrative. 



significant differences between these two sets of SSP models. 
The GALEV SSP models include the contributions of an ex- 
tensive set of nebular emission lines and gaseous continuum 
emission, which have been shown to be important in the 
first few xlO 7 yr of an SSP's evolution (Anders & Fritze-v. 
Alvensleben 2003); the Starburst99 set of SSP models does 
include nebular continuum emission, but only in a simplified 
fashion. The other main difference between both sets of SSP 
models is related to their treatment of the RSG phase, which 
is of significant importance around 10 7 yr: it is clear from 
panels (b) and (c) that using the GALEV models leaves a gap 
in the clusters' age distribution around the time that the 
RSG phase is expected to be important. This is caused by a 
combination of both the age-metallicity degeneracy and the 
sparser age resolution of the GALEV models compared to that 
of the Starburst99 SSPs. Lamers et al. (2001) have shown 
that the Geneva models of fully convective stars are not cool 
or red enough compared to the observations, particularly at 
lower metallicities (see also Massey & Olsen 2003). This may 
have significant consequences for techniques which allow the 
metallicity to be a free parameter, obviously depending on 
the age range of the respective clusters. Whitmore & Zhang 
(2002) have shown that cluster models calculated with the 
Padova tracks fit the observations better than those calcu- 
lated with the Geneva tracks. 

The relative mass distributions, shown in Fig. |7| are 
much more similar to each other than the corresponding 
age distributions when we compare the various methods and 
SSP models used. The cluster masses have all been corrected 
to the stellar IMF used by the galev SSPs (see Section 
13.2.41 . The effect of the age-metallicity degeneracy is seen 
to some extent between panels (a) and (b,c,d): as we es- 
tablished above, this degeneracy causes the cluster ages to 
be underestimated if the metallicity is overestimated (as is 
likely the case if we assume solar metallicity for the NGC 
3310 clusters), which in turn causes the cluster masses to be 
underestimated. However, this effect is minor in our NGC 
3310 cluster sample. In spite of the inherent problems of ap- 
plying KS tests to astrophysical data such as presented here, 
the results of such tests are illustrative in a comparative 
fashion: the difference in the adopted metallicity reduces the 
probability of distributions (a) and (b) to have been drawn 
from the same population to only 19.0 per cent. The agree- 
ment between distributions (b), and (c) and (d), are more 
satisfactory, with probabilities of these having been drawn 
from the same population of 93.0 and 67.3 per cent, respec- 
tively. In addition to the KS statistics, we can also compare 
the overall statistics of the distributions in Fig. |7| as shown 
in Table [3] which shows that we can reproduce the mean 
("peak") and spread of the distributions consistently and 
well within the uncertainties (represented by the a values), 
even where we adopted different metallicity distributions; if 
we simply compare the peak values of the mass distribu- 
tions obtained from the various methods, we find a spread 
among these values of <tm = A{log(M c i/Mo)) < 0.06 (where 
we have only used the peak values obtained assuming sim- 
ilar boundary conditions, i.e., for Zq\ the peak values were 
obtained from Gaussian fits to the distributions of the indi- 
vidual cluster masses). If we had simply taken the mean of 
the age distributions and done the same comparison, the re- 
sulting spread would have been at = A(log(Age/yr)) < 0.15 
(although we note that this result is unphysical, in view of 
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Figure 6. Resulting age distributions of the NGC 3310 star clus- 
ters, using a variety of methods and SSP models, as indicated in 
the individual panels. The vertical dotted lines correspond to the 
fitting boundaries of the models. Metallicities are indicated in the 
panels. 
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Figure 7. Resulting mass distributions of the NGC 3310 star 
clusters, using a variety of methods and SSP models, as indicated 
in the individual panels. Metallicities are also indicated in the 
panels. 



the significantly non-Gaussian distributions, so that a sin- 
gle "mean" value does not convey much useful information) . 
Thus, we conclude that the peaks in the relative mass distri- 
butions can be derived much more consistently than those 
in the relative age distributions. This is caused by the way 
in which the masses are determined (by scaling up either 
the entire observed SED or the T/-band flux to the mod- 
els), and the fairly narrow age range covered by the NGC 
3310 clusters (ensuring a relatively small range of cluster 
mass-to- light ratios). 

5.2 Restricted wavelength coverage, including Ha 
observations: NGC 4038/9 

While ideally one would like to have the most extensive 
wavelength coverage possible, realistically one cannot expect 
to obtain more HST coverage than by, say, four passbands 
for any given cluster system. Under this assumption, we have 
shown - both empirically (de Grijs et al. 2003a,c) and the- 
oretically (Anders et al. 2004b) - that the optical UBVI 
passbands (or their equivalents in the HST flight system) 
provide the most suitable passband combination to use as 
the basis for our broad-band SED analysis. Additional NIR 
observations would add significantly more leverage, but in 
practice such observations need to be obtained using dif- 
ferent detectors, and are thus more difficult to obtain. We 
focus therefore on data sets that can be obtained with min- 
imal observing time, while maximising the scientific output. 

In this section, we will explore the differences between 
the various methods, using UBVI coverage of 20 star clus- 
ters in the Antennae interacting galaxies (NGC 4038/9), 



selected to span a large age range in the original analysis 
where their properties were first published (see Whitmore 
et al. 1999, and references therein). The "standard" UBVI 
cluster magnitudes were obtained by converting the HST- 
flight system magnitudes using the Holtzman et al. (1995) 
conversion equations. For a subset of these clusters we have 
also obtained Ha EWs, which can - in principle - be used to 
constrain their ages more accurately and robustly, although 
strong and patchy background Ha fluxes (as observed in 
the Antennae galaxies) can render the actual contributions 
to the Ha fluxes by the clusters themselves very uncertain. 

Figure |H] shows the relative age distributions resulting 
from the application of the various methods described in Sec- 
tion^ The top three panels of Fig.|HJ (a), (b), and (c), show 
similar trends as pointed out for the same method+models 
combinations used for the NGC 3310 clusters in the pre- 
vious section. The effects of the age-metallicity degener- 
acy are somewhat less pronounced in this case, since our 
multi-dimensional AnalySED analysis (Section 13. 2. 311 indi- 
cates that the cluster metallicities in the Antennae galax- 
ies are roughly equally split between solar metallicity and 
Z ~ 2.5^0 (with only a few clusters characterised by sub- 
solar metallicities, Z < 0.01). These metallicity estimates 
are supported by high-resolution spectroscopy obtained by 
Mengel et al. (2002). Thus, by assuming solar metallicity for 
all Antennae clusters, we will have overestimated the ages 
for those clusters with supersolar metallicity, and underesti- 
mated the ages of the few subsolar metallicity clusters. This 
is reflected by the different age distributions between pan- 
els (a) vs. (b) and (c). Again, application of these methods 
leaves a clear gap in the age distributions around the age 
where the RSGs become apparent. 
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The multi-dimensional SED analysis further shows a 
weak correlation between E(S— V) and cluster age (although 
with a large scatter), from E(_B — V) < 1 mag at ~ 10 6 ' 5 yr 
to E(B - V) < 0.4 mag at ~ 10 9 ' 5 yr. 

In Fig. IHJl, we display the results from the "Sequen- 
tial O/IR" method ("Section 13.31 . Based on the availabil- 
ity of BVI photometry, this method also allows us to de- 
termine the extinction toward the sample clusters indepen- 
dently. The weak trend found by our multi-dimensional SED 
approach is also found using this method, although the 
individual extinction values are generally slightly smaller, 
E(B -V)< 0.2 mag. 

Yi et al. (2004) have shown that the (U — B) versus 
(B — V) two-colour diagram can be used to break the age- 
metallicity degeneracy for metal-poor populations. However, 
since our sample clusters are metal rich, we cannot apply this 
technique here. All of their ageing trajectories (based on the 
BC00 SSP models), whatever their metallicity, are at the 
same locus of the diagram, at least for stellar populations 
older than 100 Myr. 

While the double-peaked age distribution obtained in 
panels (a)-(c) is to some extent reproduced, the two-step 
process of the Sequential O/IR method results in a more 
evenly spread age distribution. This is partially due to the 
fact that some of the sample clusters are apparent outliers 
in the diagnostic diagrams, and as a consequence their ages 
are not well constrained. 

Finally, Figs. |H£ and f show the age distributions re- 
sulting from using the reddening-free Q-parameter analy- 
sis and the broad-band+Ha method (Sections 13.41 and 13. 51 . 
respectively. These distributions are, in very broad terms, 
consistent with those obtained using the other methods dis- 
cussed before, in the sense that they show multiple peaks 
at roughly similar ages (although they do not match in de- 
tail). The most deviant distribution is in fact that resulting 
from the "Sequential O/IR" method, which is based on a 
smaller number of photometric data points per cluster than 
the other methods. 

Figure [§] shows the corresponding mass distribution 
for our sample of Antennae clusters, based on a variety 
of methods. Because of the small-number statistics and, in 
particular, the different mass range covered by all of the 
mass distributions in Fig. |3] the simple KS statistic indi- 
cates that all distributions are different from our base line 
distributions, Figs. and b, at the > 98 per cent level 
(even if we restrict ourselves to the largest possible mass 
range in common between any two sets of mass determina- 
tions, and common metallicity assumptions). However, as 
for NGC 3310, the overall characteristics of the mass dis- 
tributions are fairly consistently reproduced, in particular 
the mean mass (see also Table similarly as for NGC 
3310, we find a spread in the mean mass among the various 
approaches of om = A(log(Af c i/M Q )) < 0.14 (once again, 
for the fits done assuming Zq only). The equivalent spread 
for the (unphysical) mean in the age distributions would be 
a t = A(log(Age/yr)} < 0.35. It appears, therefore, that - 
once more - the peaks in the relative mass distributions (as 
opposed to the detailed shapes of the distributions 4 ) can be 



4 For a comparison of the detailed shape of the underlying distri- 
bution, one would need a statistically much larger (and unbiased) 
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Figure 8. Resulting age distributions of the Antennae star clus- 
ters, using a variety of methods and SSP models, as indicated 
in the individual panels. The vertical dotted lines correspond to 
the boundaries of the models. Metallicities are indicated in the 
panels. 



obtained with a much higher degree of confidence than those 
of the relative age distributions. 

The overall characteristics of the mass distributions for 
the current sample of Antennae clusters are summarised in 
Table 2] which were obtained using the same fitting algo- 
rithms as for the fits to the NGC 3310 mass distributions. 
We note that, due to insufficient information in the Kurth et 
al. (1999) SSP models at the highest time resolution (used in 
this paper), we were unable to determine the cluster masses 
for the Q — Q method. 



6 ONE-TO-ONE COMPARISONS OF MASS 
AND AGE ESTIMATES 

Having established that the peak (or most frequent value) 
of the relative age and, in particular, mass distributions of a 
given star cluster system can be retrieved relatively robustly 
using broad-band SEDs, we will now explore to what extent 
this applies to the individual cluster properties themselves. 
It is clear from the outset that where multiple estimates for 
the cluster extinction values, and to a lesser extent also for 



sample of clusters than studied here. It should be noted that here, 
we selected clusters biased in such a way that they would cover as 
extensive a range in ages and masses as possible for these galaxies, 
based on preliminary analysis (see Section |5J. For smaller sam- 
ples such as ours, the uncertainties in the individual age and mass 
estimates (which are on the order of the histogram bin sizes; cf. 
de Grijs et al. 2003a, c, Anders et al. 2004b) start to affect the re- 
sulting distribution non-negligibly and in ways that cannot easily 
be quantified. 
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Figure 9. Resulting mass distributions of the Antennae star clus- 
ters, using a variety of methods and SSP models, as indicated in 
the individual panels. Metallicities are also indicated in the pan- 
els. 



Table 4. Characteristics of the overall relative mass distributions 
of our Antennae cluster sample 
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their metallicities, exist, these estimates vary significantly. 
However, despite this being the case from an observational 
point of view, the effects on the SED of varying the ex- 
tinction and metallicity, if their uncertainties are not more 
than a few tenths in magnitude or one step in metallicity, 
respectively, are small, and can thus still result in reason- 
ably secure age and mass estimates for all but the youngest 
clusters (cf. Anders et al. 2004b). In other words, the age 
and mass estimates are relatively robust to variations in the 
extinction and metallicity estimates. Therefore, we will only 
discuss the one-to-one comparisons of the ages and masses of 
the individual clusters in either of our samples. We also note 
that the youngest SEDs (< 10 7 yr) can be affected quite sig- 
nificantly (and not in any systematic fashion) by even small 
changes in extinction and/or metallicity, however (see An- 
ders et al. 2004b). 

Figure tTTjl shows the one-to-one comparisons of the NGC 
3310 cluster ages and masses, adopting solar metallicity, and 
using the AnalySED approach as our basis for the compari- 
son. It is immediately clear that the methods using the same 



set of SSPs, Figs. llOh and c for the ages and masses, respec- 
tively, result in highly reproducible age and mass estimates 
(well within the 1-a error bars). With few exceptions, the in- 
dividual ages and masses match very well within the model 
uncertainties. 

In the most extreme case considered for the NGC 3310 
star cluster system, namely by using a different method 
(AnalySED vs. 3DEF) and a different set of SSP models 
(galev vs. Starburst99; shown in Figs.llUb and d), the indi- 
vidual age and mass estimates still match up well within the 
model uncertainties, although with a slightly larger scatter. 
The "3DEF,Starburst99" approach results in slightly higher 
masses than the AnalySED approach, although this is only 
a ~ la deviation. 

Figure 1111 shows a similar set of comparisons for the 
Antennae clusters, now using the entire range of methods 
and models at our disposal. At first sight, we notice three 
characteristics when comparing these panels to the panels of 
Fig. 1101 First, the uncertainties in the ages and masses esti- 
mated by the AnalySED approach (see Anders et al. 2004b 
for a full description and justification) are significantly larger 
than those resulting from most other methods used, except 
for the 3DEF approach. Secondly, the uncertainties are much 
greater than those for the individual age and mass estimates 
obtained for the NGC 3310 clusters. This reflects the differ- 
ence in wavelength coverage of the broad-band SEDs for 
both sets of clusters: with the smaller wavelength coverage 
of the Antennae clusters, their best-fitting ages and masses 
are not as well constrained as for the more extensively cov- 
ered NGC 3310 clusters. It therefore appears that the un- 
certainties estimated using the Sequential O/IR, the Q-Q 
and the BB+Ha methods are too small in view of the ill- 
defined broad-band SEDs over the fairly small wavelength 
range covered. We should point out, however, that the uncer- 
tainties resulting from the BB+Ha method (i.e., for clusters 
with measured Ha EWs) are likely to be smaller than those 
from the broad-band-only methods, since this method uses 
Ha EWs to further constrain the cluster ages. This method 
is therefore most accurate in a very narrow age range (i.e., 
more accurate than methods based on broad-band photom- 
etry alone), around ~ 10 Myr. Finally, it appears that the 
mass estimates of the individual clusters are better matched 
than their age estimates. This, in turn, provides the more 
robust relative mass distributions discussed in the previous 
Sections. 

Thus, we believe that the main differences among the 
resulting age and mass estimates for individual clusters are 
caused by the difference in the wavelength range covered, 
but - as we will show in Section|7|- the uncertainties intro- 
duced by transforming the iTST-flight system magnitudes to 
"standard" ground-based UBVI photometry may affect the 
robustness of the results to a similar, if not greater, degree. 
While the individual cluster ages and masses based on the 
"ideal" U BVI (or equivalent) coverage advocated in de Grijs 
et al. (2003c), and Anders et al. (2004b) may be subject 
to significant uncertainties, the relative age and mass dis- 
tributions are much more consistently established, so that 
statistical analyses of large-scale cluster systems based on 
multi-passband broad-band imaging do have the potential 
to provide robust scientific insights in their host galaxies' 
formation, evolution, and star-formation histories. 

Finally, in Table [S] we compare the age estimates for the 
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Figure 10. One-to-one comparisons between the age (left) and 
mass (right) estimates of the NGC 3310 clusters, for solar mctal- 
licity, obtained from the various methods and SSP models, using 
our AnalySED approach as the basis for comparison. The dotted 
lines are the loci of equality. 
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Figure 11. One-to-one comparisons between the age (left) and 
mass (right) estimates of the Antennae clusters, for solar mctal- 
licity, obtained from the various methods and SSP models, using 
our AnalySED approach as the basis for comparison. The dotted 
lines are the loci of equality. 



individual clusters in our Antennae cluster sample with Ha 
EW > l.OA, for the full set of SSP analysis methods used. 
We have indicated the clearly discrepant age estimates in 
italic font. The final two columns of this table include the 
mean age and its standard deviation based on the individ- 
ual age determinations; for those clusters with clearly dis- 
crepant values, we also give these numbers excluding these 
discrepant determinations. In addition, we have also in- 
cluded those cluster for which we determined ages of log( 
Age/yr ) < 6.6 but for which no strong Ha EW measure- 
ments were provided. 

One can see immediately that the correspondence be- 
tween the age estimates for a given cluster among the dif- 
ferent analysis approaches is close. We note that the age es- 
timates obtained using the sequential O/IR and AnalySED 
approaches, with metallicity as a free parameter are offset 
with respect to most of the other methods. This simply re- 
flects the age-metallicity degeneracy for these young ages. 
Nevertheless, the fact that we obtain very similar age esti- 
mates using a variety of independent approaches that may 
or may not include Ha luminosities as an additional con- 
traint is reassuring. In principle, this validates the differ- 
ent approaches used here: nearly all clusters that should be 
young based on the presence of strong Ha emission are in- 
deed young (with between and 2 outliers for all methods), 
and the ones predicted to be young based on these methods 
all have Ha (except for, again, between and 2 outliers) . The 
outliers are characterised by the largest error bars, which is 
an additional argument in support of the robustness of the 
variety of models employed here. Unfortunately, there are 
no young extragalactic star clusters available in the current 
literature for which independent age estimates have been ob- 
tained via either spectroscopy or detailed analysis of their 
resolved colour-magnitude diagrams; the use of Ha EWs, as 
done here, is therefore the closest we can get to an indepen- 
dent validation of our approach. 

As an example of the uncertainties inherent to the use 
of broad-band (and Ha) fluxes to obtain the individual clus- 
ter ages, we direct the reader's attention to cluster G2-14, 
for which we found the most discrepant age estimates among 
the variety of approaches used. This cluster shows strong Ha 
emission, and must therefore be young. Nevertheless, two 
of the approaches employed assigned this cluster an age of 
greater than 1 Gyr. Upon close inspection, those discrepant 
estimates originated from the two approaches that essen- 
tially leave the metallicity as a free parameter. This hints at 
an origin related to the age-metallicity degeneracy. However, 
we also point out that the error bars assigned by the Anal- 
ySED approach arc among the largest in our cluster sample. 
Statistically, the AnalySED age estimate (with metallicity 
as a free parameter) is therefore also consistent with a young 
age. 

On the other hand, all of the different approaches em- 
ployed here, for instance, estimate the age of cluster G2-17 
to be in the range 6.6 < log(Age/yr) < 6.9 (see the Ap- 
pendix for details), despite not having strong Ha emission. 
In spite of this, the uniformity of this result among the dif- 
ferent methods is encouraging; it implies a high degree of 
reproducibility. 
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Table 5. Age comparison for the individual clusters in our Antennae cluster sample with EWh q > 1.0A. 
Values in italics are clearly discrepant values. 







Q-Q 


BB+Hq 


AnalyoHilJ; ^ tree 




.lySED; Zq 




O /IR 


log( Agc/yr ) 






log( Ag 


;c/yr ) 








Mean 




G2-04 


2.9 


6.5 


6.7 


7.1 




6.9 


6.8 


7.2 


6.87 


0.24 


G2-05 


2.5 


7.5 


6.8 


7.1 




6.6 


6.8 


7.6 


6.98 


0.29 




















(6.83 


0.18) 6 


G2-07 


1.3 


6.5 


6.5 


7.1 




6.6 


6.8 


7.1 


6.77 


0.26 


G2-08 


2.9 


6.5 


6.8 


7.2 




6.6 


6.8 


7.0 


6.82 


0.23 


G2-10 


1.4 


6.6 


6.8 


6.9 




6.9 


6.8 


7.4 


6.90 


0.24 


G2-12 


2.5 


6.6 


6.8 


7.3 




6.6 


6.8 


6.8 


6.82 


0.23 


G2-14 


3.7 


6.5 


6.1 


6.9 




6.6 


9.0 


9.3 


7.40 


1.26 




















(6.53 


0.29) 6 


G2-15 


3.5 


6.5 


6.5 


7.1 




6.6 


6.6 


6.8 


6.68 


0.21 


Others 


with 




















log( Af 


;c/yr ) < 6.6, 


G2-17 


G2-19 


G2-17 a 




G2-09, 




G2-17 a 






but no 


Ha 










G2-19 











NOTES: a Here, we used a limit of log( Agc/yr ) < 6.9 because of the offset introduced by the age-metallicity degeneracy; After exclusion of the values in 
italic font. 



7 FILTER SYSTEM CONVERSIONS 

Largely owing to the unavailability of SSP models computed 
for the HST-flight system magnitudes, many of the early 
studies based on HST imaging observations of extragalactic 
star cluster systems used the equations given by Holtzman et 
al. (1995) to convert HST STMAG magnitudes to the "stan- 
dard" Johnson- Cousins system. Despite recent updates of 
many of the leading SSP models, which now include theo- 
retical magnitudes in the STMAG system, many workers in 
the field, including ourselves in this paper (see our photome- 
try of the Antennae clusters) , continue to use the Holtzman 
et al. (1995) conversions. 

It is clear, however, that any conversion based on 
generic spectral properties of a given stellar population will 
introduce biases and additional uncertainties that could, in 
principle, be avoided by retaining one's photometry in the 
filter system used for the observations. In this section, we 
explore the extent of these additional uncertainties by com- 
paring our model fitting results for the full NGC 3310 cluster 
sample, presented in de Grijs et al. (2003c), based on both 
the original HST/WFPC2 photometry and on the trans- 
formed magnitudes in the Landolt (KPNO) UBVRI system 
used by Holtzman et al. (1995). 

In order to do so, we folded the GALEV SSP models 
through these filter transmission curves, which were kindly 
made available by Jon Holtzman (priv. comm.). Despite be- 
ing a standard system, the appropriate filter transmission 
curves have not been published in their entirety. We trans- 
formed our F336W, F439W, F606W and F814W magnitudes 
using Eqs. (6)-(9) below (for the WF3 chip of the WFPC2 
camera), using an iterative approach; following Holtzman 
et al.'s (1995) recommendations, we did not transform the 
F300W magnitudes, but used the HST flight system mag- 
nitude for this part of the broad-band cluster SEDs. The 
conversions from the F336W, F439W, and F814W niters 
to their UBI counterparts are based on the observational 
transformations from the WFPC2 flight system (Holtzman 
et al.'s Table 7), while the transformation of the F606W 
magnitudes to the V filter relies on the conversion of the 
synthetic WFPC2 system to the standard V band (their Ta- 
ble 10). 



U = -2.5 x log(counts s _1 ) - 0.240 x (U - V) + 

0.048 x(U- V) 2 + 18.764 + 2.5 x log(2.003) (6) 
B = -2.5 x log(counts s^ 1 ) + 0.003 x (B — V) + 

-0.088 X (B — Vf + 20.070 + 2.5 x log(2.003) (7) 
V = -2.5 x log(counts s _1 ) + 0.254 x (V - I) + 

+0.012 x (V - I) 2 + 22.093 + 2.5 x log(2.003) (8) 

(for (V — I) < 2.0) 
I = -2.5 x log(counts s _1 ) - 0.062 x (V - I) + 

0.025 x (V - I) 2 + 20.839 + 2.5 x log(2.003) (9) 

As before, we will first discuss the characteristics (i.e., 
predominantly the mean and spread) of the age and mass 
distributions of the entire cluster population, and conclude 
with a one-to-one comparison of the results obtained for 
the individual clusters. For the analysis presented in this 
section, we have adopted the AnalySED multi-dimensional 
modelling approach, using the GALEV SSP models. We at- 
tempted to obtain the cluster ages and masses under three 
sets of assumptions: (i) unrestricted fits, i.e., we left all of 
the cluster ages, masses, metallicities and extinction values 
as free parameters; (ii) as for (i), but assuming solar metal- 
licity for all clusters; (Hi) as for (ii), but now also assuming 
a generic (arbitrarily low) extinction value for each cluster 
of E(B -V) = 0.1 mag. 

In Fig. 1121 we present the results for the relative age 
distributions, in the left-hand column using the HST-flight 
system magnitudes as our basis, and in the right-hand col- 
umn using the transformed F300W- UB VI photometry. From 
top to bottom, the fits become more and more restricted, fol- 
lowing the assumptions laid out above. Despite this being 
the same sample as analysed in de Grijs et al. (2003a,c), the 
resulting age distribution in panel (a) is different from that 
published previously. The main reason for this difference is 
that the fits discussed in this section do not include the NIR 
passbands, so that we are less sensitive to metallicity varia- 
tions (Anders et al. 2004b). In other words, the strong peak 
seen at log(Age/yr) ~ 6.6 is caused by the age-metallicity 
effect, and also by the fact that the youngest age in our 
models corresponds to log(Age/yr) = 6.6 (so that younger 
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clusters will automatically be assigned the minimum age in 
the models). This applies to the strong peaks seen at this 
age in all of the panels of Fig. 1121 We note that the overall 
metallicity of NGC 3310 is known to be significantly sub- 
solar (see Section I5.1H . but we will use the assumption of 
solar metallicity in this section to emphasize a number of 
technical concerns relevant for similar studies in this field. 

Apart from this obvious signature of the age-metallicity 
degeneracy, the overall relative age distribution, i.e., the 
mean value and its spread, based on the HST flight sys- 
tem photometry is retrieved relatively robustly under the 
various assumptions employed (see Fig. 1121 left-hand pan- 
els); the range of ages found for the NGC 3310 clusters is 
consistent with independently determined age estimates for 
the star-forming regions in this galaxy (see de Grijs et al. 
2003a,c for a detailed discussion) . However, if we now exam- 
ine the age distributions based on the transformed F300W- 
UBVI magnitudes, we see that (i) our results are severely 
affected by the age-metallicity degeneracy, and (ii) that the 
"transformed ages" do not correspond to even remotely sim- 
ilar ages as obtained from the HST-flight magnitudes. Only 
by severely restricting our model fits (Fig. I12b 2'l are we able 
to retrieve a similar age distribution as we obtained from 
the direct match of our SSP model grid to the iTST-flight 
system magnitudes. This provides, therefore, a very strong 
argument against using photometry based on even robust 
filter conversions; the effect of such transformations is the 
unavoidable introduction of biases and additional uncertain- 
ties, which thus makes the various fitting routines less reli- 
able and robust. 

If we now consider the resulting mass distributions, 
shown in Fig. 1131 we see - somewhat to our surprise - that 
to first order these can be reproduced relatively robustly, 
based on either set of filter transmission curves. This is ow- 
ing to the fact that our mass estimates are predominantly 
determined by a global scaling between the entire observed 
SED and the most appropriate model SED, rather than on 
the exact shape of the SED. 

To highlight the robustness with which we can retrieve 
the global characteristics of the mass distribution, and also 
to address the effects caused by the age-metallicity degen- 
eracy, we show the relative mass distributions for only the 
clusters with ages, log(Age/yr) > 6.8 as the cross-hatched 
histograms, where relevant. The global characteristics, in- 
cluding the mean and spread (a), of all mass distributions 
shown in Fig. 1131 are listed in Table |S| We note that while 
the global characteristics of the mass distributions are repro- 
duced robustly, the details of the distributions differ among 
the panels of Fig. 1131 To illustrate this, we applied KS tests 
to the relevant data sets, of which the results are summarised 
in Table □ 

Finally, we compare the individual cluster age and mass 
estimates obtained from both sets of filter transmission 
curves, and all three sets of assumptions in Fig. 1141 The indi- 
vidual panels in this figure reflect the discussion above: the 
individual age estimates are severely discrepant, while the 
mass estimates are relatively robust from one cluster to an- 
other. We note the existence of two "sequences" in Figs. ll4H 
and e. These are indeed caused by the age-metallicity degen- 
eracy discussed above; the upper ( "left-hand") sequence con- 
sists predominantly of clusters that are found in the strong 
peak at our minimum age limit. 




67896789 



Log( Age/yr ) log( Age/yr ) 

Figure 12. Comparison of the relative age distribution of the 
NGC 3310 clusters, based on the HST flight system photometry 
(left-hand column, panels al, bl, and cl) and on the converted 
UBVI magnitudes (right-hand column, panels a2, b2, and c2). 
Panels al and a2 are based on model fits in which all of the 
cluster ages, masses, metallicities and extinction values were left 
as free parameters; in panels bl and b2 we restricted the fits to 
solar metallicity, and in panels cl and c2 we also adopted a generic 
extinction of E(B — V) = 0.1 mag. 



Table 7. Detailed comparison of the mass distributions in Fig. 
1131 Probabilities refer to the chances that both the "HST" and the 
"UBVI" samples were drawn from the same distribution, based 
on KS statistics. 



Selection 


Restrictions 


Probability 


all 


unrestricted 


0.663 


log( Age/yr) > 6.8 


unrestricted 


0.379 


all 


z e 


0.250 


all 


Z Q , E(S - V) = 0.1 


0.047 


log(Ago/yr) > 6.8 


Z e , B(B - V) = 0.1 


0.015 



Thus, we conclude from this analysis that the overall 
characteristics of a cluster system's relative mass distribu- 
tion, and to some extent the individual cluster mass esti- 
mates as well, can be reproduced fairly robustly under a 
variety of relevant fitting assumptions and conversions be- 
tween filter systems. However, in order to derive more robust 
age estimates, one should ideally retain one's photometry in 
the filter system used for the observations. This is in order 
to introduce as few biases and additional uncertainties as 
possible. We emphasize that we have followed this route in 
previous publications in which we applied these results. 
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Table 6. Characteristics of the overall relative mass distributions of the NGC 3310 cluster sample, as discussed w.r.t. Fig. 1131 respectively. 



System 


Selection 




Restrictions 




log(Ma 
Mean 


5 s/M ) 
a 


HST 


all 




unrestricted 




5.03 


0.41 


HST 


log(Age/yr) > 6.8 




unrestricted 




5.17 


0.33 


HST 


all 




z 




4.92 


0.44 


HST 


log(Age/yr) > 6.8 




z 




5.14 


0.29 


HST 


all 




E(B - V) = 


0.1 


4.87 


0.46 


HST 


log(Agc/yr) > 6.8 


z e 


E(S - V) = 


0.1 


5.20 


0.24 


UBVI 


all 




unrestricted 




5.11 


0.40 


UBVI 


log(Age/yr) > 6.8 




unrestricted 




5.18 


0.39 


UBVI 


all 




z& 




4.87 


0.36 


UBVI 


all 


z e 


E(B - V") = 


0.1 


5.08 


0.50 


UBVI 


log(Age/yr) > 6.8 


z @ 
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0.1 


5.38 


0.26 
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Figure 13. Comparison of the relative mass distribution of the 
NGC 3310 clusters, based on the HST flight system photometry 
(left-hand column, panels al, bl, and cl) and on the converted 
UBVI magnitudes (right-hand column, panels a2, b2, and c2). 
Panels al and a2 are based on model fits in which all of the 
cluster ages, masses, metallicities and extinction values were left 
as free parameters; in panels bl and b2 we restricted the fits to 
solar metallicity, and in panels cl and c2 we also adopted a generic 
extinction of E(B— V) = 0.1 mag. The cross-hatched distributions 
contain clusters characterised by log( Age/yr ) > 6.8, to avoid the 
effects of the age-metallicity degeneracy (see text). 

8 SUMMARY AND CONCLUSIONS 

The increasing availability of high-resolution HST imaging 
observations across a wide wavelength range (from the mid- 
UV to the NIR) has revolutionised studies of extragalactic 
star cluster populations. 

The age, mass, metallicity and extinction of unresolved 
extragalactic star clusters can be derived from the spectral 
energy distributions measured in broad-band photometric 
systems, by comparing them with cluster evolution models. 
This method is applied in the literature based on observa- 



Figure 14. One-to-one comparisons between the age (left) and 
mass (right) estimates of the NGC 3310 clusters, based on the 
HST flight system photometry vs. the converted UBVI magni- 
tudes. The dotted lines are the loci of equality. 

tions with different sets of photometric filters and different 
sets of cluster evolution models. 

In this paper, we investigated the accuracy of these de- 
terminations by comparing the parameters derived from dif- 
ferent filter sets and different sets of cluster models. To this 
purpose, we used the UBVI and HST photometry of 20 
clusters in the Antennae galaxies and 17 clusters in NGC 
3310, analysed in the different ways that are found in the 
literature. We assess the systematic uncertainties in age and 
mass determinations, and to a lesser extent also in extinc- 
tion and metallicity determinations. We compare the results 
with those of the extensively tested and well- validated Anal- 
ySED approach (see Sections 13.71 and HI . The results, which 
are summarised below, give us a handle on the systematic 
uncertainties one needs to contend with when comparing 
results obtained by different groups, each using a different 
modelling technique. 
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We first examined the model parameters used by the 
various groups that contributed to the final results. We con- 
cluded that the differences among the various SSP models, 
specifically among the various Bruzual & Chariot model in- 
carnations, the Starburst99 and the GALEV SSP models, are 
random and do not bias the results in a systematic fash- 
ion. Similarly, the variety of (foreground) extinction curves 
used by the various groups show minimal differences over the 
wavelength range considered in this paper (from the mid-UV 
to the NIR), and any differences are deemed unimportant in 
view of the photometric uncertainties of similar magnitude. 
Thus, we conclude that any significant differences among 
the resulting parameters are due to the details of the var- 
ious methods used, rather than to the models themselves 
(see also Yi 2003). 

The methods used in this study are: 

(i) the multi-dimensional SED analysis (Anders et al. 
2004b) and (Bik et al. 2003); 

(ii) the broad-band fluxes with Ha method (Whitmore & 
Zhang 2002); 

(iii) the reddening-free parameter method (Q — Q) (Whit- 
more et al. 1999), and 

(iv) the optical/NIR sequential analysis method (Par- 
mentier et al. 2003). 

The free parameters in the matching of cluster mod- 
els to an observed photometric cluster SED are, mass, age, 
metallicity and extinction. 

The effects of (foreground) extinction and metallicity 
variations are small, within a given cluster sample. However, 
on an individual cluster basis, it is very difficult to estimate 
a best-fitting metallicity, despite major efforts and recent 
improvements made in the modelling techniques (e.g., de 
Grijs et al. 2003a,b, Anders et al. 2004b). This agrees with 
the conclusion by Bastian et al. (2004) that they could not 
determine the metallicity of individual clusters in the in- 
teracting spiral galaxy M51, but that the derived age and 
mass distributions are only slightly affected compared to the 
solar-metallicity case. Thus, the effect of metallicity varia- 
tions is minimal in view of the uncertainties, and conclusions 
regarding cluster age and mass distributions appear robust. 

For this reason, and also in order to compare our results 
to previously published results, we adopted solar metallici- 
ties for most of the modelling done in this paper. We find 
that the ages of a cluster, measured using different meth- 
ods, can differ drastically, but that the "age distribution" 
of a given young (< 10 9 yr) cluster system derived using 
the different methods for a fixed metallicity shows the same 
main features. We also find that the differences in the mass 
distributions, derived using different approaches, are much 
smaller than those in the age distributions. Thus, the mass 
distributions can be obtained with a higher degree of confi- 
dence than the age distributions. We determine accuracies 
for our age and mass distributions, based on the simplify- 
ing assumption that these properties are roughly following 
a Gaussian distribution, of o t = A(log(Age/yr)} < 0.35 and 
om = A(log(M c i/M©)) < 0.14, respectively. While this as- 
sumption may not hold in general, the differences resulting 
from its blanket application to all results presented in this 
paper are indicative of the degree of confidence we can at- 
tribute to our fits. The very small spread in the mean mass 
shows that, if there is a peak in the mass distribution, the 



retrieval of its absolute value is relatively insensitive to the 
approach taken. 

Since the actual age distributions of our cluster sam- 
ples are distinctly not unimodal, the use of a "mean" value 
and "width" for the comparison of methods provides a 
straightforward, yet not a physically interesting criterion. 
The absolute age distribution is rather sensitive to the 
adopted method, and one should therefore be very cautious 
when comparing such values among different studies. As we 
showed previously (e.g., de Grijs et al. 2001, 2003a,b,c, An- 
ders et al. 2004a,b), the relative age distributions (i.e., the 
presence or absence of peaks) can be retrieved with a very 
high level of confidence, however. 

It is important to keep in mind that all of the techniques 
employed in this paper are affected by various selection ef- 
fects and artifacts that require careful consideration when 
comparing results. For instance, the peak at young ages in 
the age distribution of Fig. 1131 is caused by the use of a 
model of this age as the youngest time-step, while the bi- 
modality in a number of the resulting age distributions is 
due to differences in the treatment of RSG stars (see also 
Whitmore & Zhang 2002). Nevertheless, the close correspon- 
dence among the resulting cluster parameters based on the 
variety of methods used is reassuring in the context of the 
robustness of our parameter determination techniques. 

We note that as extensive a wavelength coverage as pos- 
sible is required to obtain robust age and mass estimates 
for the individual objects, with reasonable uncertainties. We 
also show that the conversion of the HST photometry to the 
"standard" Johnson and Cousins UBVI photometry intro- 
duces an extra uncertainty in the determination of the ages 
and masses of the clusters, because the calibration of this 
conversion is based on stars rather than clusters. 
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Table A9. Results for the Antennae clusters from the Q-Q ap- 
proach, using the Kurth et al. (1999) SSP models and adopting 
solar metallicity for all clusters (cf. Fig. 151:1 . 



Table A10. Results for the Antennae clusters from the BB+Ha 
approach, using the BCOO SSP models and adopting solar metal- 
licity for all clusters (cf. Figs.|HJ and |9*H . 



ID 


A v (mag) 


l°g( Agc/yr ) 


ID 


E(S - V) (mag) 


log( Agc/yr ) 


Mass (xlO 5 M ) 


G2-01 


0.0 


9 . 3 


G2-01 


().()() 


8.56 


0.96 


G2-02 


0.0 


9.3 


G2-02 


0.00 


8.61 


1.03 


G2-03 


0.8 


8.0 


G2-03 


0.38 


8.01 


0.49 


G2-04 


1.3 


6.5 


G2-04 


0.34 


6.79 


0.21 


G2-05 


1.9 


7.5 


G2-05 


0.76 


6.78 


0.45 


G2-06 


0.0 


9.3 


G2-06 


0.00 


8.41 


0.77 


G2-07 


0.6 


6.5 


G2-07 


0.42 


6.46 


0.003 


G2-08 


2.4 


6.5 


G2-08 


0.60 


6.76 


0.25 


G2-09 


1.6 


8.0 


G2-09 


0.78 


6.68 


0.66 


G2-10 


1.6 


6.6 


G2-10 


0.42 


6.84 


0.48 


G2-11 


2.8 


7.7 


G2-11 


1.10 


6.64 


0.26 


G2-12 


0.8 


6.6 


G2-12 


0.08 


6.78 


0.64 


G2-13 


2.8 


7.9 


G2-13 


0.08 


9.38 


1.29 


G2-14 


4.9 


6.5 


G2-14 


1.20 


6.10 


1.08 


G2-15 


1.7 


6.5 


G2-15 


0.44 


6.46 


1.47 


G2-16 


0.4 


7.7 


G2-16 


0.44 


6.68 


0.007 


G2-17 


1.2 


6.6 


G2-17 


0.30 


6.86 


0.66 


G2-18 


0.0 


9.3 


G2-18 


0.00 


8.61 


0.74 


G2-19 


2.3 


7.2 


G2-19 


0.96 


6.58 


0.20 


G2-20 


2.2 


9.0 


G2-20 


0.24 


9.16 


1.92 



APPENDIX A: RESULTS FROM THE 
INDIVIDUAL APPROACHES 

In this Appendix, we present the best-fitting results for 
the individual clusters as obtained by applying the different 
methods discussed in this paper to the cluster photometry 
presented in Section [5] The following tables are arranged 
following the order of the age and mass histograms in Figs. 

nan 
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Table Al. Results for the NGC 3310 clusters from the multi-dimensional SED fits, using the AnalySED approach with cluster ages, 
masses, metallicities and extinction values as free parameters (cf. Figs.|§k and 17k .) . 
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z 






E(S 




V) ( 


mag 
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M 


;c (x 10 7 


yr) 
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3 (XlO 5 


Mq) 
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Mm. 


Best 


Max. 


Mm. 


Best 


Max. 


Min. 


Best 


Max. 


Gl-01 
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.0004 





0005 


0.17 





.25 


0. 


32 


0.84 


1 


.20 


1 


.56 


1.88 


2 


.69 


3.50 


Gl-02 





.0004 





.0004 


0. 


0004 


0.00 





.00 


0. 


10 


1.20 


1 


.60 


1. 


(id 


1.95 


1 


.95 


2.23 


Gl-03 





.0040 





0080 





0080 


0.00 





.00 


0. 


15 


1.60 


3 


.60 


4 


80 


0.55 





.73 


0.94 


Gl-04 





.0200 





0200 





0200 


0.00 





.05 


0. 


05 


1.20 


2 


.40 


2 


.80 


0.27 





.66 


0.77 


Gl-05 





.0004 





.0040 


0. 


0080 


0.00 





.10 


0. 


30 


2.80 


6 


.00 


12 


40 


0.78 


1 


.07 


1.83 


Gl-06 





.0003 
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0005 


0.07 





.10 


0. 


13 


0.84 


1 


.20 


1 


.56 


1.40 


2 


.00 


2.60 


Gl-07 





.0004 





0004 





0200 


0.20 





.20 


0. 


25 


0.40 





.80 





.80 


2.67 


5 


.62 
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00 
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0500 
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00 


2.80 
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.80 
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.80 


0.40 





.58 


0.70 
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.0004 
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2.') 
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2.80 


Gl-11 





.0040 





0040 





0200 


0.00 





.15 
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.00 
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0.61 





.99 


1.45 


Gl-12 





.0040 





.0040 


0. 


0080 


0.00 





.05 


0. 


05 
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05 


4.40 


5 


.20 


12 


00 


0.74 





.70 


1.10 
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.0040 
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0.00 





.15 


0. 


35 


1.20 


4 


.80 


23. 


20 


0.49 





.89 
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0. 


0200 
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80 
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0200 


0.00 





.55 


0. 


65 


0.80 
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.60 


60. 


40 


0.17 


2 


11 


3.10 


Gl-17 





.0200 





0500 





0500 


0.00 





00 


0. 


00 


5.60 


5 


60 


13. 


20 


1.00 
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.01 


1.60 



Table A2. Results for the NGC 3310 clusters from the multi-dimensional SED fits, using the AnalySED approach with cluster ages, 
masses, and extinction values as free parameters, and adopting solar metallicity for all clusters (cf. Figs. 03 andlTb). 



ID 




B(S 




V) ( 


mag 


) 




Age (XlO 7 


y) 






Mass 


(XlO 5 


m q) 


Min. 


Best 


Max. 


Min. 


Best 


Max. 


Min. 


Best 


Max. 


Gl-01 





.17 





.25 


0. 


32 





28 





.40 





.52 





.47 





67 


0.87 


Gl-02 


0. 


07 





.10 





13 


0. 


28 





.40 





.52 





.39 


0. 


55 


0.72 


Gl-03 


0. 


00 





.00 





00 


5. 


20 


5 


.00 


5. 




1 


.09 


1 


17 


1.17 


Gl-04 





00 





.00 





05 





80 





.80 


2 


00 





.27 





27 


0.67 


Gl-05 


0. 


.05 





.05 





10 


4. 


00 


6 


.00 


(i 


0(1 





.75 


1 


03 


1.17 


Gl-06 


0. 


10 





.15 





19 


0. 


28 





.40 





.52 





.41 


0. 


58 


0.76 


Gl-07 


0. 


17 





.25 





32 





28 





.40 





.52 


1 


.87 


2. 


67 


3.47 


Gl-08 





.35 





.40 





40 


0. 


40 





.40 





40 





.13 


0. 


10 


0.16 


Gl-09 


0. 


00 





.00 





00 


6. 


40 


(i 


.40 


6 


.80 





.66 


0. 


68 


0.68 


Gl-10 





.21 





.30 


0. 


39 





28 





.40 





.52 





.57 


0. 


82 


1.06 


Gl-11 


0. 


00 





.00 





00 


2. 


80 


3 


.20 


4 


40 





.81 





88 


1.09 


Gl-12 


0. 


17 





.25 





32 





28 





.40 





.52 





.92 


1 


31 


1.70 


Gl-13 





00 





.00 


0. 


00 


7 


2(1 


7 


.00 


7 


(i() 





.88 





88 


0.88 


Gl-14 


0. 


.10 





.15 





15 


3. 


60 


5 


AM 


6 


40 





.79 


1 


14 


1.20 


Gl-15 





00 





.05 





05 


4. 


40 


5 


.00 


6. 


0(1 





.82 


1 


09 


1.14 


Gl-16 


0. 


00 





.00 





25 


5. 


20 


30 


.00 


48 


.80 


1 


.08 


1 


56 


2.44 


Gl-17 


0. 


00 





.00 


0. 


00 


7 


20 


7 


.00 


8 


40 


1 


.27 


1 


28 


1.32 



Table A3. Results for the NGC 3310 clusters from the 3DEF approach, using the GALEV SSP models and adopting solar metallicity for 
all clusters (cf. Figs.[Sk and 17k). 



ID E(B - V) (mag) log( Agc/yr ) log( Mass/M Q ) 





Min. 


Best 


Max. 


Min. 


Best 


Max. 


Min. 


Best 


Max. 


Gl-01 


0.22 





26 


0. 


32 


6.60 


(i 


(i() 


6 


(id 


4.65 


4 


7(1 


4.78 


Gl-02 


0.06 


0. 


12 


0. 


Hi 


6.60 


6 


(id 


6 


(id 


4.62 


4 


69 


4.74 


Gl-03 


0.00 


0. 


00 


0. 


04 


7.51 


7 


.78 


7 


.81 


4.81 


5 


(12 


5.07 


Gl-04 


0.00 


0. 


04 


0. 


Hi 


7.08 


7 


.08 


7 


.56 


4.41 


4 


45 


5.02 


Gl-05 


0.00 





10 


0. 


18 


6.90 


7 


.78 


7. 


.86 


4.d6 


4 


.98 


5.05 


Gl-06 


0.06 


0. 


12 


0. 


18 


6.60 


(i 


(id 


6 


(id 


4.54 


4 


63 


4.71 


Gl-07 


0.24 


0. 


28 


0. 


32 


6.60 


6 


(id 


6 


(id 


5.29 


5 


34 


5.38 


Gl-08 


0.00 


0. 


00 


0. 


10 


6.90 


7. 


.81 


7 


.86 


3.7d 


4 


(id 


4.72 


Gl-09 


0.00 





00 


0. 


10 


7.38 


7 


.83 


7 


.88 


4.50 


4 


.82 


4.92 


Gl-10 


0.26 


0. 


32 


0. 


36 


6.60 


(i 


(id 


6 


(id 


4.74 


4 


.81 


4.86 


Gl-11 


0.00 


0. 


06 


0. 


12 


6.90 


7. 


(id 


7 


.86 


4.24 


5 


.03 


5.25 


Gl-12 


0.22 





28 


0. 


34 


6.60 


6 


(id 


6 


(id 


4.97 


5 


04 


5.10 


Gl-13 


0.00 


0. 


00 


0. 


14 


7.45 


7. 


.88 


8 


.13 


4.71 


4 


.95 


5.08 


Gl-14 


0.00 


0. 


18 


0. 


30 


6.90 


7 


.56 


8 


.43 


4.11 


4 


.84 


5.22 


Gl-15 


0.00 





08 


0. 


18 


6.90 


7. 


.81 


7 


.88 


4.06 


5 


.05 


5.16 


Gl-16 


0.00 


0. 


30 


0. 


48 


6.90 


7 


.81 


8 


.95 


4.13 


5 


(19 


5.51 


Gl-17 


0.00 





00 


0. 


12 


7.86 


8 


09 


8 


25 


5.07 


5 


17 


5.28 
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Table A4. Results for the NGC 3310 clusters from the 3DEF approach, using the Starburst99 SSP models and adopting solar metallicity 
for all clusters (cf. Figs. 01 andlTlll. 



ID E(B - V) (mag) log( Agc/yr ) log( Mass/M ) 





Min. 


Best 


Max. 


Min. 


Best 


Max. 


Min. 


Best 


Max. 


Gl-01 


0.00 


0. 


36 





42 


6.20 


6 


41 


7.28 


4.41 


4 


.54 


4.65 


Gl-02 


0.06 


0. 


16 





2(1 


6.56 


6 


56 


6.61 


4.34 


4 


41 


4.46 


Gl-03 


0.00 


0. 


04 





16 


6.82 


7 


.19 


7.82 


3.70 


4 


27 


4.77 


Gl-04 


0.00 


0. 


00 





10 


6.85 


7. 


.12 


7.15 


3.75 


4 


.16 


4.32 


Gl-05 


0.00 


0. 


08 





26 


6.82 


7 


90 


8.34 


3.54 


4 


.78 


5.08 


Gl-06 


0.08 


0. 


16 





26 


6.20 


(i 


56 


6.61 


4.26 


4 


.35 


4.52 


Gl-07 


0.28 


0. 


32 





36 


6.41 


6. 


.56 


6.56 


5.01 


5 


06 


5.14 


Gl-08 


0.00 


0. 


00 





46 


6.56 


7 


.98 


8.37 


3.26 


4 


45 


4.75 


Gl-09 


0.00 


0. 


18 





26 


6.78 


7. 


.23 


7.67 


3.51 


4 


.27 


4.46 


Gl-10 


0.00 


0. 


34 





42 


6.20 


6 


49 


7.28 


4.45 


4 


50 


4.66 


Gl-11 


0.00 


0. 


04 





.58 


6.20 


7. 


.15 


7.98 


3.78 


4 


.28 


4.97 


Gl-12 


0.00 


0. 


00 





.38 


5.00 


7 


.24 


7.32 


4.66 


4 


90 


4.95 


Gl-13 


0.02 





12 





.32 


6.78 


7 


44 


7.71 


3.62 


4 


49 


4.75 


Gl-14 


0.00 


0. 


18 





70 


5.78 


7 


45 


8.16 


3.52 


4 


.45 


5.02 


Gl-15 


0.00 


0. 


08 





.58 


6.49 


7 


72 


8.30 


3.62 


4 


70 


5.10 


Gl-16 


0.00 


0. 


40 





.92 


5.00 


7. 


.28 


8.77 


3.54 


4 


50 


5.41 


Gl-17 


0.14 


0. 


20 





.38 


6.78 


7. 


.44 


7.58 


4.15 


4 


08 


4.81 



Table A5. Results for the Antennae clusters from the multi-dimensional SED fits, using the AnalySED approach with cluster ages, 
masses, metallicities and extinction values as free parameters (cf. Figs.|SJi and l9*h~) . 



E(B - V) (mag) Age ( X 10° yr) Mass ( X 10° Mq) 

[in. Best Max. Min. Best Max. Min. Best Max 



G2-01 





0004 





0500 





0500 





00 





20 





40 





.56 


1 


.76 


10.50 





.55 


1 


63 


2.47 


G2-02 


0. 


0004 





0080 





.0500 





00 





.15 





45 





4S 


3. 


.96 


11.70 





.54 


1. 


.84 


2.54 


G2-03 





0004 





0500 
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Table A6. Results for the Antennae clusters from the multi-dimensional SED fits, using the AnalySED approach with cluster ages, 
masses, and extinction values as free parameters, and adopting solar metallicity for all clusters (cf. Figs. 03 and 19b - ). 



E(B - V) (mag) Ago ( X 10 8 yr) Mass ( X 10 J Mq ) 
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Table A7. Results for the Antennae clusters from the 3DEF approach, using the Starburst99 SSP models and adopting solar metallicity 
for all clusters (cf. Figs.^; and lot). 



ID E(B - V) (mag) log( Agc/yr ) log( Mass/M ) 
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Table A8. Results for the Antennae clusters from the "Sequential O/IR" approach, using the BCOO SSP models and metallicities ranging 
from 0.4 to 2.5Zq, as indicated in the column headings; the M/Ly ratios are given for solar metallicity (cf. Figs.[5Ji and|5Jl). 
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4.89 


G2-07 





00 
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4.98 


G2-09 


0. 
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